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Abstract 



We present the general theory of clean, two-dimensional, quantum Heisenberg 
antiferromagnets which are close to the zero-temperature quantum transition 
between ground states with and without long-range Neel order. While some 
of our discussion is more general, the bulk of our theory will be restricted 
to antiferromagnets in which the Neel order is described by a 3-vector order 
parameter. For Neel-ordered states, 'nearly-critical' means that the ground 
state spin-stiffness, p s , satisfies p s <C J, where J is the nearest-neighbor 
exchange constant, while 'nearly-critical' quantum-disordered ground states 
have a energy-gap, A, towards excitations with spin-1, which satisfies A -C J. 
The allowed temperatures, T, are also smaller than J, but no restrictions are 
placed on the values of ksT/p s or ksT/A. Under these circumstances, we 
show that the wavevector/frequency-dependent uniform and staggered spin 
susceptibilities, and the specific heat, are completely universal functions of 
just three thermodynamic parameters. On the ordered side, these three pa- 
rameters are p s , the T = spin- wave velocity c, and the ground state stag- 
gered moment Nq; previous works have noted the universal dependence of the 
susceptibilities on these three parameters only in the more restricted regime of 
ksT <C p s . On the disordered side the three thermodynamic parameters are 
A, c, and the spin-1 quasiparticle residue A. Explicit results for the universal 
scaling functions are obtained by a 1/N expansion on the O(N) quantum 
non-linear sigma model, and by Monte Carlo simulations. These calculations 
lead to a variety of testable predictions for neutron scattering, NMR, and 
magnetization measurements. Our results are in good agreement with a num- 
ber of numerical simulations and experiments on undoped and lightly-doped 
La2-sSrsCu04. 
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I. INTRODUCTION 



The subject of two-dimensional quantum antiferromagnetism has witnessed a remark- 
able revival in recent yearsEHl. This is largely due to the intense interest in understanding 
the properties of the CuO<i layers in the high temperature superconductors. However, the 
experimental motivation is not limited to these cuprate compounds; recent investigational 
have refocussed interest on a number of other layered insulating compounds which are rather 
well described as Heisenberg antiferromagnets at low temperatures^!. 

The bulk of the existing theoretical work on two dimensional antiferromagnets can be 
divided into two broad classes: 

First, there are the studies of the low temperature properties of antiferromagnets with 
well-established long-range Neel order in their ground stateS'i@0. Definitive results have 
been obtained for these systems by Chakravarty eia/.9il. They showed that the long- 
wavelength, low energy properties were well described by a mapping to a classical two- 
dimensional Heisenberg magnet. All effects of quantum fluctuations could be absorbed into 
almost innocuous renormalizations of the coupling constants. Good agreement with neutron 
scattering experiments on La^CuO^ was obtained. 

Second, there have been numerous investigations on spin-fluid, or quantum disordered, 
ground statesliUll. These are states in which quantum-fluctuations have removed all vestiges 
of the Neel order. An entirely new physical picture is necessary for visualizing the ground 
states: it is phrased most often in terms of resonating singlet valence-bonds between pairs 
of quantum spins0 Many interesting questions on the presence of alternative symmetry 
breaking in the ground state have been addressed. The nature of the excitations above the 
ground state is also of some interest. The two experimentally distinguishable possibilities 
are (i) the low- lying states correspond to those associated with weakly-interacting spin- 
1/2 quanta ('spinons'), or {%%) the spinons are confined in pairs, leading, to integer-spin 
excitations with infinite lifetimes. The reader is referred to a recent reviewa where many of 
these questions are addressed in greater detail. 

In this paper, we present a detailed theory of quantum antiferromagnets which fall in 
between the above two classes. These nearly-critical antiferromagnets are neither strongly 
Neel ordered nor fully quantum disordered in the ground state. The ordered Neel moment, 
if it exists, is much smaller than the ordering moment of the corresponding classical antifer- 
romagnet. If the ordering moment is absent, the resulting quantum-disordered ground state 
has an energy-gap towards excitations with non-zero spin which is much smaller than all 
microscopic energy scales in the problem. Such nearly critical antiferromagnets were con- 
sidered briefly by Chakravarty et. aZ.0; they identified three different regimes of behavior 
(See Fig. [I]) which we now outline. At short length and energy scales the spin correlations in 
these antiferromagnets are essentially critical. The spins fluctuate strongly between ordered 
and non-ordered configurations. At very low temperatures, the quantum-fluctuating system 
only makes up its mind at a fairly large scale, and crosses over to behavior characteristic of 
either a Neel-ordered ground state (this is the 'renormalized-classical' region of Fig. [XJ) or a 
quantum-disordered ground state (the 'quantum-disordered' region of Fig. ^ see also Fig. 0). 
At larger T there is another, very interesting possibility: the critical quantum fluctuations 
may be quenched by thermal effects before the system has had a chance to undergo the above 
crossover to one of its two ground states (the 'quantum-critical' region of Figs. and||). The 
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system then does not display the properties of either Neel-order or quantum-disorder at any 
length scale; instead it crosses over from critical spin fluctuations to a thermally-induced, 
quantum-relaxational regime which will be described for the first time in this paper. 

The above three crossover occur at a large length-scale in nearly-critical antiferromagnets, 
suggesting that their properties should be universal. The bulk of this paper is devoted to 
making this statement more precise. We will indeed find that the long-wavelength, low- 
energy uniform and staggered spin susceptibilities are completely characterized by universal 
scaling functions. The contribution of the spins to the specific heat will also be found to be 
similarly universal. The only required inputs are three thermodynamic parameters, which 
will be described more precisely below. Some of the results reported here have been discussed 
briefly in recent reports!!!. Some other recent work has also discussed quantum criticality 
near related magnetic phase transitionsllEI. 

Our motivation in examining nearly-critical antiferromagnets comes primarily from nu- 
merous recent experiments on undoped and weakly doped La2-sSr$CuO±. 

We consider first undoped Tc^CmC^. It is by now well established that La^CuO^ is de- 
scribed extremely well as a spin-1/2 square lattice Heisenberg antiferromagnet with nearest- 
neighbor interactions^. Almost all theoretical studies@i!c! of this system have focussed 
primarily on its properties at very low temperatures, where a description in terms of classi- 
cal fluctuations of the Neel order parameter is appropriate; this range of temperatures was 
referred to as the 'renormalized-classical' regionH There is good accord between theoryS'ElIl 
and experiments!! - !! at these low temperatures (T): the correlation length £(T) increases 
exponentially with falling T, the equal-time structure factor, S(k), at zero momentum be- 
haves as S(0) oc T 2 £ 2 , and the 63 Cu spin-lattice relaxation rate 1/Ti decreases rapidly with 
the increasing T, all of which is in good agreement with the renormalized classical theory. 
However, recent experiments^ have shown that at intermediate temperatures (T > 0.4J 
where J is the nearest-neighbor exchange constant), 1/Ti becomes nearly independent of T. 
The crossover to this behavior occurs at a T which is small compared to J, so one can hope 
that a low-energy theory of the expected crossover to the quantum-critical region might 
still be appropriate; the possibility of such a crossover was already noted earlieiO. We will 
show in this paper that this new behavior is well described quantitatively by a theory of the 
quantum-critical spin fluctuations. The correlation length £(T) is also expected to display 
a crossover at these temperatures^; unfortunately, there are no experimental data for £ for 
T > 0.4J. However, strong support for our interpretation comes from the experimental 
and numericalil~@ measurements of the uniform static spin susceptibility x?(T). We will 
show that its temperature-dependent slope above T ~ 0.35 J is in excellent accord with the 
predictions of the quantum-critical theory, and in clear disagreement (by a factor of three) 
with the result deduced from a renormalized-classical theory. Taken together, we will argue 
that the above results imply the following: the square lattice spin-1/2 Heisenberg antiferro- 
magnet with nearest-neighbor exchange has long-range Neel order in its ground state, and is 
well described by a renormalized classical theory at low temperatures; however it is appar- 
ently close enough to a quantum phase transition to a quantum disordered phase to display 
quantum-critical spin fluctuations over an appreciable range of intermediate temperatures. 

Consider next weakly doped La^-^Sr^CuO^. Nonzero doping (5 > 0) has two important 
consequences: (i) the bare value of the spin-stiffness p s takes a smaller value, pushing the 
antiferromagnet closer to the quantum phase transition, and (ii) the antiferromagnet is per- 
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turbed by the random potential of the dopant ions. Randomness is a relevant perturbation 
at the T = quantum phase transitioned and must be included in any theory of the low 
T properties. Dynamic neutron scattering measurementsllS near 5 = 0.04 show that the 
spin spectrum can be collapsed in a scaling plot in which the measurement frequency u is 
scaled by T; a related dependence of the susceptibility on uj and T was also discussed in 
the phenomenological theory of the marginal Fermi liquidlll. We have argued that this uj /T 
scaling is in fact a rather general property of quantum-critical spin fluctuations, even in the 
presence of randomness and dopingcl. Furthermore, the pre-factor of the scaling function at 
low uj and T shows clear evidence of the importance of randomnessS However many other 
experimentsiiB on doped La 2 -sSrsCuOi have been performed at relatively large tempera- 
tures. For these T we assume that the antiferromagnet is insensitive to the weak randomness 
and is still in the vicinity of the pure fixed point. The primary effect of non-zero doping will 
then be to reduce the value of p s - an immediate consequence is that the NMR relaxation 
1/Ti should be nearly T independent (the quantum-critical behavior) over a T range which 
increases with doping. This is indeed what is observed. We will discuss this and other 
comparisons with experiments on the doped cuprates in more detail later in Section | VII . 



With a single adjustable parameter (the doping dependent p s ) we will find good agreement 
between our theory and the experiments. 

We turn now to a more complete description of our results. Consider the antiferromagnet 
described by the following Hamiltonian 

TC = 22JijS i -S j , (1.1) 

i<j 

where i,j extend over the sites of two-dimensional lattice, Sj are on-site spin operators 
acting on states with spin S on the site i, and the Jy are exchange constants which fall off 
rapidly with the separation between i and j. The Jij are assumed to be invariant under the 
translation symmetry of the underlying lattice. The energy scale J will be used to denote 
the largest of the The are predominantly antiferromagnetic, so that the classical 
ground state has no average uniform magnetization. We will be interested primarily in 
antiferromagnets which undergo a zero-temperature quantum phase transition from a Neel- 
ordered to a quantum disordered state as the ratios of the are varied, and the strength of 
the quantum fluctuation increases. Let us represent this strength by an all-purpose coupling 
constant g; the system is assumed to be Neel-ordered for g smaller than a critical coupling g c 
and quantum-disordered for g > g c . We will assume further that the critical ground state at 
g = g c is described by a continuum field theory of excitations which propagate with a non- 
singular spin-wave velocity c - a number of explicit mean-field solutions of such transitions 
have this property^E, although other possibilities have also been discussedSSS. It should 
be possible to extend our results to antiferromagnets which have different velocities for 
different spin-wave polarizations00, but we will not consider this complication here. In 
Section |VII| we will argue that the following scaling results also apply unchanged to the 
corresponding quantum phase transitions in lightly-doped antiferromagnets. 

The nature of the classical ordering helps us identify the proper continuum fields nec- 
essary for a hydrodynamic theory of the quantum phase transition. The first of these is of 
course the Neel order parameter. We will restrict the analysis in this paper to antiferromag- 
nets with an ordinary vector order parameter. This type of order parameter is associated 
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with ground states with collinear spin-ordering i.e. the on-site spin condensates are either 
parallel or anti-parallel to each other. More complicated order parameters can arise in sys- 
tems with stronger frustration e.g. the triangular and kagome lattices which have co-planar 
spins and matrix order parameters^lHB; we will not discuss these complications here. Re- 
turning to the vector order-parameter case, we assume that the condensate is oriented in the 
rkz direction, and define a continuum quantum field n(r), which will be used as the hydro- 
dynamic order-parameter variable, as: n(rj) = Sj on sites where the condensate points up, 
and n z {vi) = —Sf] n^C^i) = on sites where the condensate points down. The staggered 
magnetization N is then 

N = (n z (r)) T=0 . (1.2) 
Upon approaching the critical point at T = 0, this staggered magnetization will vanish 

N ~(g c -g) /3 , (1.3) 

where (3 is a universal critical exponent. We have iVo = for g > g c . Furthermore, 
(n z (r)) = at all finite T because it is not possible to break a continuous non-abelian 
symmetry in a two-dimensional system. At the critical point, g = g c , equal-time n(r) 
correlations will decay with an anomalous power lawll 

WO • n(0)> ~ (1.4) 

where D = 3 is the dimension of space-time. 

The second important hydrodynamic variable is the magnetization density quantum field 
M(r) 

M(r<) = ^S,, (1.5) 
cr 

where a 2 is the volume per spin, and g/iB is the gyromagnetic ratio of each spin. Although 
the ordering has no net magnetization, magnetization fluctuations decay slowly due to the 
conservation law on the total magnetization. 

Finally, the Hamiltonian TC is itself associated with a conserved total energy. The contri- 
bution of the spins to the specific heat per unit volume, Cy, is the appropriate experimental 
observable, sensitive to this hydrodynamic quantity. 

The hydrodynamic properties of the order parameter and magnetization fluctuations can 
be determined from the following two retarded response functions: 



Xs{k,uj)6 i>m = —^J ^ r J Q dt 



([n,(r,t) ) n m (0,0)])e-^ r -^ ) 

^ p pco 

Xu(k,uj)5e m = J d 2 r J dt 

([Mt(r, t), M m (0, 0)])e" i(k ' r ^* ) , (1.6) 
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where all fields have now acquired a Heisenberg-picture time (t) dependence, the indices £, m 
extend over the three spin directions x, y, z, and the average is with respect to a thermal 
Gibbs ensemble at a temperature T. These correlation function are the dynamic staggered 
and uniform spin-susceptibilities respectively; their values predict the result of essentially 
all the experiments that have been performed on antiferromagnets. An important exception 
is the Raman-scattering cross-sectionE! - we will not discuss its properties here. 

We now present the scaling forms satisfied by Xs, Xu, arid Cy in the vicinity of the 
quantum phase transition at g — g c . The temperature, T, is taken to be non-zero, but must 
satisfy 

k B T < J. (1.7) 

A non-zero T implies the absence of a spin-condensate, and the response functions are 
therefore rotationally invariant. It is useful to describe separately the scaling properties 
of magnets with and without Neel order in their ground state. This will be followed by a 
discussion of the relationship between the two cases. 



A. Neel ordered ground state: 

The scaling properties should clearly depend upon a variable which measures the distance 
of the ground state from criticality. The most convenient choice is the ground state spin- 
stiffnessEl p s . Its value can be easily determined by experiments and by various numerical 
analyses on model Hamiltonians, with no arbitrary overall scale factors. In two-dimensions, 
p s has the dimensions of energy, and the requirement that the magnet is not too far from 
criticality is 

Ps < J. (1.8) 
Upon approaching g c , p s obeys Josephson scaling 

Ps~{g c -g) {D - 2) \ (1.9) 

where v is the usual correlation length exponent. We can now state one of the central results 
of this paper: For g < g c , and under the conditions on T and p s noted above, the values of 
Xsi Xu and Cy satisfy the following scaling forms 

„ s A n 2 / he V ( Nk B T\ v ( hck huo Nk B T\ 

^'"Hid ^ t$i H^t'^t'^7J' (L10b) 

3C(3), (k B T\ 2 (Nk B T\ 

where iV is the number of components of the order-parameter, £ is the Reimann zeta function, 
and <3>i s , $i„ and are completely universal, dimensionless, functions ($i u and $i s are 
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complex while is real) denned such that they remain finite as T/p s — > oo; this will 
also be true for other scaling functions introduced below. For simplicity, we have explicitly 
specialized to antiferromagnets with spacetime dimension D = 2 + 1, although analogous 
results for general D are not difficult to write down. Particularly striking is the absence 
of any non-universal scale factors (in either the arguments or the prefactors of the scaling 
functions) in all scaling forms. Everything is fully determined by the values of p s , c and Nq 
and there is no further dependence on lattice scale physics. The universal dependence of the 
spin susceptibilities on p s , N and c was implicit in the analysis of Chakravarty et. a/.IH (see 
also the recent work of Hasenfratz and Niedermayeii!) for the low T regime T C p s ; our 
results are however valid for all values of T/p,. Also Castro Neto aad EradidrH have re.eer.tly 
discussed closely related scaling forms for Cy near general quantum phase transitions in 2 + 1 
dimensions. The coefficient of \&i in ( |1.10c| ) has been chosen to be the specific heat of a 
single gapless bose degree of freedom with dispersion u = ck in 2 dimensions. The number 
^>i(T — > 0) is thus a measure of the effective number of such modes in the ground state. 
This number is given by \l/i(0) for the ordered Neel phase, and by \Pi(oo) for the quantum 
critical state at g = g c . 

We note further that all scaling forms continue to be valid even at g — g c : the prefactor 
of $ ls in ( P-.10a| ) remains non-singular at g — g c because of the results (|1.3|), ( |1.9| ) and the 
exponent identity!*! 

2(3 = (D -2 + r])u. (1.11) 

The arguments of the scaling functions will occur frequently in this paper. We therefore 
introduce the dimensionless variables 

— Tick huj 

which represent momentum and frequency measured in units of a scale set by the absolute 
temperature T. The third argument 

Nk B T 

Xl = —^-, 1.13 
2np s 

determines whether the antiferromagnet is better described at the longest distances as a 
quantum-critical or a renormalized-classical model (see Fig. [I]). The factor of iV in the 
definition of X\ is to facilitate the large N limit in which p s ~ N; the variable X\ will therefore 
remain of order unity. The factor of l/(27r) is purely for future notational convenience. For 
large x±, the energy scale k B T is the largest energy which first cuts-off the critical spin 
fluctuations, and the system never fully realize that its coupling g is in fact different from 
g c and that the ground state is ordered: the spin-fluctuations are quantum-critical at the 
shortest scales, and are eventually quenched in a universal way by the temperature. For small 
X\ the antiferromagnet is in the renormalized-classical region. There is a large intermediate 
scale over which the antiferromagnet behaves as if it has long-range Neel order; eventually, 
strong two-dimensional classical thermal fluctuations of the order-parameter destroy the 
Neel order. 
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B. Quantum-disordered ground state: 



We now consider the case g > g c . We will assume that the quantum-disordered state 
has a gap towards all excitations. This has certainly been satisfied by all explicit large N 
constructions of such states in frustrated antiferromagnets in the vicinity of the transition 
to long-range Neel orderEHil. We will assume further that there are no deconfined spin-1/2 
excitations above theground state: this is expected to be true in systems with a collinear 
Neel order parametert3. Antiferromagnets on the triangular or kagome lattices with coplanar 
spin correlations are expected to possess deconfined spin-1/2 spinon excitations0@ - their 
critical properties will therefore not be described by the present theory. For the case of 
confined spinons which is under consideration here, the lowest excitation with non-zero spin 
will carry spin 1. Further, at T = 0, this excitation should have an infinite lifetime at small 
enough k. The distance from criticality is conveniently specified by the gap, A, to this spin 
1 excitation. The equal-time order parameter correlation function will decay exponentially 
on a scale £ which is inversely proportional to A. Therefore A will vanish as 

A~(g-g c y (1.14) 

upon approaching criticality. One is not too far from criticality provided 

A < J. (1.15) 

Further, we need an observable which sets the scale for order-parameter fluctuations. On 
the ordered side this was done by N . A convenient choice on the disordered side is to use 
an amplitude of the local, on-site, dynamic susceptibility, xl- This susceptibility is defined 
by 

i f°° 

XL{u)$em = -T dt([Su(t),S irn (0)]) 
n Jo 

r d 2 k 

~ Sim J -^Xs(k,u). (1.16) 



In principle, Xu also contributes to xl, but when ( |1.15| ) is satisfied its contribution is sub 



dominant to that from Xs, and can therefore be neglected. For g > g Cl it can be shown that 
at T = 0, xl has the following imaginary part for small u> close enough to the threshold A: 



lm XL (oo)\ T=0 = j sgn(u;)0(%| - A), (1.17) 



where 9 is the unit step function, and A/ 4 is an amplitude with the dimensions of inverse- 
energy. We will show later that the discontinuity iA/A in the local dynamic susceptibility 
is precisely a quarter of the quasiparticle residue A of the low-lying spin-1 excitation. As g 
approaches g c , A vanishes as 

A-ig-g^. (1.18) 

We have now assembled all the variables necessary for obtaining the scaling forms for Xs 
and Xu for g > g c . The relations analogous to ( |1.10a| ), ( |1.10b| ) and ( |1 . 10c| ) are 
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M = [ ^ )kBT ^ { ^^^\, (1 . 19b) 




Or A. (*-£)' (1.19c) 

where $2s, 3*2^5 and \l/2 are completely universal functions, which are finite in the limit 
T/A — > 00. The physical response functions are again completely determined by three 
thermodynamic parameters: A, c, and A, with no further sensitivity to lattice scale physics. 
As in Sec |I^, all scaling forms continue to be valid even at g — g c : the prefactor of $2s in 
( |1.19a| ) remains non-singular at g — g c because of the scaling results (|1.14Q , ( J1.18| j . 
We also introduce, for future convenience, the variable 

x 2 = (1 .2o) 

which determines whether the antiferromagnet is in the quantum-critical or quantum dis- 
ordered regions (see Fig [I]). For large x 2 , the temperature T predominates the small zero- 
temperature gap, A, and the system may as well be at g = g c . For small x 2 , the ground-state 
gap A quenches the spin fluctuations, putting the system in the quantum-disordered region. 
The thermodynamics is well described in terms of a dilute gas of activated excitations. 



C. Critical point: 

We have obtained above two separate universal scaling forms at the critical coupling 
g = g c , but T finite, by taking the limits g S g c and g \ g c . It follows therefore that the 
two results must be simply related: 

$ 2s (k,u, xi = 00) = Z Q § ls (k,L}, x 2 = 00), 

§2u(k,Uj,Xi = OO) = $ lu (k,Lj,X 2 = OO), 

* 2 (oo) = * 1 (oo), (1.21) 

where Zq is a universal number. Recall that the universal functions have been chosen to 
have a a finite limit as x\$, — > 00. 

Note that there is no rescaling factor for the uniform susceptibility and the specific heat. 
This is because their overall scale is universal and was not set by some thermodynamic 
observable, as was the case for the staggered susceptibility. This universality in scale is 
related to the fact that M and the energy are conserved quantities: this will discussed 
further in Sec y. For the staggered susceptibility, we have performed a 1/N expansion of 
the scale factor on antiferromagnets with an iV-component order-parameter and found 

0.229191243 

Zq = 1 ^ • (1.22) 
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Actually it is not just the values, but the entire asymptotic expansions of the universal 
functions which have matching conditions at x\$ = oo. As the antiferromagnet has no phase 
transition at finite temperature, all its measurable properties should be smooth functions of 
the bare coupling constant g — g c provided T^O. This fact, combined with x\ ~ (g c ~g) v , 
x 2 ~ (g — g c )~" ', can be used to easily deduce the constraints on the asymptotic expansions 
of $u, $m, ^1, ®2s, $2«, ^2 about x 1;2 = oo. 



D. Experimental Observables 

The main purpose of the rest of the paper is to describe the universal functions $i n , 
$2s; ^2u, and ^2 as completely as possible. An large amount of information is contained 
in them; in particular, as we shall see later, in a suitable limit they contain the complete 
static and dynamic scaling functions of Chakravarty et. a/.0 and Tyc et. alM. A large 
number of experimentally testable quantities can be obtained from these functions; now we 
highlight some of the most important by endowing them with their own scaling functions. 

We begin with the measurements related to Xu- 



1. Static, uniform spin susceptibility and spin diffusivity 

The conservation of M makes the small k and u dependence of \ u rather simple. In the 
hydrodynamic limit, ut£ 1, where T£ is a typical lifetime of excitations, the magnetization 
fluctuations must obey a diffusion equation; we find for g < g c that 

— Di (xi )k 2 — 

$i M (fc,o7,Xi) = Qi(zi) _ V ; =2; fc,57«l, (1.23) 

— iu + Di(x\)k 

and an analogous expression for g > g c with 1-^2. The functions Qi(xi), ^2(^2), -Di(^i), 
D 2 (x 2 ) are all universal. From ( |1.10b| ), ( |1.12| ), ( |1.19b| ) we see that they are related to the 



static, uniform spin susceptibility xj, and spin diffusion constant D s . We have for g < g c 

Ds(T) = ^A(^i), (1.24) 
k B T 

and analogous expressions for g > g c with 1 —>■ 2. 



The Wilson ratio is defined by 



2. Wilson Ratio 
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Its properties are therefore easily obtainable from our scaling results for \u an d Cy. It 



follows from (|1.10cj) , (|1.19c|) and ( 1.24 ) that W is a completely universal function of x\ {X2) 
for g <g c {g> g c ). We have for g < g c 



W 



7T Qi(xi) 

303) ^i( Xl y 



1.26) 



and analogous expression for g > g c with 1 — > 2. 
We turn next to experiments sensitive to \s- 



3. Structure factor 
The equal-time spin structure factor, S(k), is related to Xs(k,uj) by 

S{k)S em = I C / 2 r(n £ (r,0)n m (0,0))e- ik r 
roc du 5, 



£t i-X ^ b '- (L27) 



From (|1.10a|) , we deduce that for g < g c it satisfies the scaling form 



N 2 (fir) 2 

^ = —^^£1(^1), (1-28) 

where the universal function Si is given by 

/oo djjj \ _ 
— ^=1m$ u (k,W,xi). (1.29) 
-co 7T 1 — e 

Similarly for g > g c , we can relate S(k) to a universal function S 2 (fc, £2) with the prefactor 
Nq/p s replaced by A and the subscript 1 — ► 2. 



^. Antiferromagnetic correlation length 

We define the correlation length £ from the long-distance, e~ r ^, decay of the equal time 
n - n correlation function. This correlation function will have such an exponential decay 
for all g provided T 7^ (the actual asymptotic form also has powers of r as a prefactor). 
Equivalently, one can define £ as where in is the location of the pole of S(k) closest to 
the real k axis. The scaling function for £ for g < g c is 

r 1 = (i-so) 

For g > g c we have an identical form with 1 — > 2. Since the correlation lengths must match 
at g — g c , we clearly have Xi(oo) = X 2 (oo). The universal linear T dependence of at 
9 — 9c ( x i — 00) w as noted by Chakravarty et. alM. 
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5. Local susceptibility 



From its definition ( J1.16 ), and the scaling form ( Jl.lOa ), we can deduce the following for 
9 < 9c- 

N 2 

haxdu) = -^x , l\uJ\ r 'F 1 (uJ,x 1 ), (1.31) 

Ps 

where the universal function F\ is 

i^foari) = ^ / P^I™$is(k,uJ,x 1 ). (1.32) 
of 1 J 4tt z 

On general grounds we expect Ijxixl(w) ~ u for small uo and T non-zero; this implies that 

Fi ~ sgn(cJ) IcZ;! 1 " 77 for small UJ. (1.33) 

In principle, the real part of xl also has a singular piece which satisfies a scaling form 
analogous to Q1.31]) ; however the momentum integral in ( |1.32|) is divergent at the upper 



cutoff (provided rj > 0). Thus Re%L has a leading contribution which is non- universal and 
dominated by lattice scale physics. There is no such problem for IthXl however - in this 
case the momentum integral sums over intermediate states which are on-shell and only long- 
wavelengths contribute. Finally we note that a similar scaling form for Imx^ for g > g c can 
be obtained by replacing N$/ p s by A and the substitution 1 — > 2. 



6. NMR relaxation rate 

We consider the relaxation of nuclear spins coupled to electronic spins of the antifer- 
romagnet (e.g. Cu spins in La2CuO^. We assume that the relaxation is dominated by 
contributions near the antiferromagnetic ordering wavevector. After suitably accounting for 
lattice-scale form factors and integrating out high-energy lattice excitations, a coupling A n 
between the nuclear spins and the antiferromagnetic order parameter n can be obtained. 
The typical frequencies in NMR experiments are much smaller than the temperature and 
the relaxation rate, 1/T\ of the nuclear spins is given by 

h 1 ^ 2 ^!^^ (L34) 

From (|1.31| ) we deduce the following result for \jT\ for g < g c 



i 942/^2 

^ = ^L x v Rl ( Xl ), (1.35) 
1 1 rip s 



where i?i(xi) is a universal function given by 



iM^) = Hm (1-36) 



12 



As before, the scaling form for 1/T\ for g > g c involves replacing N^/p s by A and replacing 
1 — > 2 on the right-hand-side. 

We will discuss the general form of our results for the scaling functions for the different 
regions of the phase diagram in turn (Fig. []]). We will consider first the quantum-critical 
region (xi 1, X2 1), followed by the renormalized-classical (x\ <C 1) and the quantum- 
disordered (22 <§C 1) regions. Precise numerical results will not be presented here: the reader 
is referred to Sections pTI[ - [VI| for precise results obtained in the 1/N expansion of the O(N) 
non- linear sigma model. 



E. Quantum-critical region (x\ > 1 or 12 > 1) 

At short length/time scales the spin fluctuations in any nearly-critical antiferromagnet 
should be indistinguishable from those at the critical point. The special property of the 
quantum-critical region is that the deviations from criticality at longer length/time scales 
arise primarily from the presence of a finite T. The fact that the ground state of the system is 
not exactly at the critical point is never terribly important, and the system does not display 
behavior characteristic of either ground state at any length/energy scale. The critical spin 
fluctuations are instead quenched in a universal way by thermal relaxational effects. 

It should therefore be evident that there are two distinct types of spin fluctuations in 
wavevector /frequency space (Fig. |^). With either Tick or Tiu significantly larger than ksT, 
the spin-dynamics is that of the critical 2+1 dimensional field theory of the critical point 
g = g c . Otherwise, damping from thermally-excited, critical spin-waves produces a regime 
of quantum-relaxational dynamics. 

The crossover from the 2+1 critical to quantum-relaxational behavior is clearly evident 
in the dynamic staggered susceptibility (See Fig. |3|). In the 2+1 dimensional critical region 
(k 3> 1 or uj ^> 1) we have 

— An — 

$ ls (k,oJ, 00) = I fc»lora;»l, (1.37) 

(k -uj 2 ) 1 ^/ 2 

where Aq is a universal number, and the exponent r\ = 8/(37r 2 iV) at order 1/N but is 
knownS to order 1/iV 3 For N = 3, precision Monte Carlo simulations^ place the value 
of 77 around rj « 0.028; this extremely small, but positive, value of 77 will have important 
consequences for experiments. Note that in this regime Im$ ls is non-zero only for uJ > k 
where one obtains a broadband spectrum of critical spin-waves. The small value of 77 implies 
however that the damping is small and the spectrum is almost a delta-function. In the 
quantum-relaxational regime (k -C 1 and UJ -C 1), there is strong damping due to thermally- 
excited, critical spin-waves and excitations are not well-defined. However the spectrum of 
overdamped spin-waves remains universal and is shown in Fig. |3| 

The same crossover is also present in the structure factor. Its universal scaling function 
00) = ZQ 1 E 2 (k, 00) has the 2+1 critical form Ei(fc, 00) ~ k +v at large h, and a 
Lorentzian form at small k. 

The various static observables have a value set by the absolute temperature in a universal 
way. Their T dependence can be deduced easily from the scaling forms with the knowledge 
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that all scaling functions were chosen to have a finite limit as x± t 2 — > oo. The first corrections 
awa4y from x li2 = oo are given by0 

U{x 1>2 ) =U oc + U 1>2 /x\(Z + ... x 1)2 > 1 (1.38) 

where U represents any of the scaling functions X, Q, \1/ for the correlation length, uniform 
static susceptibility, and specific heat respectively. The form of the subleading term above 
follows from the requirement that the physics at finite temperature is a smooth function of 
the bare coupling g — g c . Chakravarty et. a/.0 have noted the X\ = oo result for the case of 
the correlation length. 



F. Renormalized classical region {x\ <C 1) 

We now describe our results at low temperatures on the ordered side, g < g c , Xi_^ 1- 
Extensive results on a closely related regime have been obtained by Chakravarty et. a/.LLHl and 
Tyc et. The relationship between our and their results is discussed below. We will find 
that, in the appropriate limit, our scaling functions reduce exactly to theirs. 

We begin with a qualitative discussion of the nature of the spin correlations in the 
renormalized-classical region. The spin-fluctuations now fall naturally into three different 
regimes in wavevector- frequency space (see Fig. 0). At the largest k,uJ we have 2+1 dimen- 
sional critical spin fluctuations which are essentially identical to those in the quantum critical 
region and are therefore described by a staggered spin susceptibility $i s similar to that in 
QL37| ). Upon moving to longer distances/times, the first crossover occurs at length (time) 
scales of order £j (£j/c) to a 'Goldstone' regime where the spin dynamics is well described by 
rotationally averaged spin-wave fluctuations about a Neel ordered ground state. The scale 
£j, controlling the critical to Goldstone crossover, is the Josephson correlation length^, and 
determines the vicinity of the ground state of the antiferromagnet to the quantum phase 
transition. Near g c , £j diverges as 

Sj~(9c-9)-"- (1-39) 

The second crossover in the renormalized classical region (Fig. |2|) occurs at the length scale 
£, which is the actual correlation length. At this scale, strong, classical, two-dimensional, 
thermal fluctuations of locally Neel ordered regions destroy the long range Neel order, so that 
at scales larger than £, the antiferromagnet again appears disordered, with all equal-time 
spin correlations decaying exponentially in space. The scale £ is roughly given by 

where we have omitted pre-exponential power- law factors of X\. For small xi, £ is clearly 
much larger than £j; the three regimes in Fig. ^| are therefore well separated. 

Our explicit results for the scaling functions will be restricted to the vicinity of the second 
crossover described above: thus they are valid when x\ <C 1, k£j 1, and uj^j/c 1. In 
this regime, all of our scaling functions, which in general depend upon three arguments k, uJ, 
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and xi, collapse into reduced scaling function depending only on two arguments measuring 
momentum and frequency, and oo^/c, Further, the reduced scaling functions turn out 
to be exactly those obtained by Chakravarty et. alM Tyc et. alHl and Hasenfratz et. 
a/.EiB. At first sight this result may seem a bit surprising. Our results were obtained 
for antiferromagnets with a small stiffness, i.e. p s <C J, while in other workBSBS, no 
restrictions was placed on the value of p s . The equivalence to order 1/N between the two 
theories is a consequence of the fact that the low T results in the renormalized classical 
region contain no corrections of order p s /J to order 1/N. In a recent analysis Hasenfratz 
et. a/.Hil. have suggested that such corrections are in fact absent at low T at all orders in 
1/N. 

We show the nature of the collapse in the scaling functions explicitly for the structure 
factor. Chakravarty et. aZ.IHI proposed the scaling form 

S(k) = 5(0) /(*£), (1.41) 

where / is a universal function. We find that / is related to the scaling function Si by 

f(y) = A/ lim o {( Xi rV^Xl( Xl ) 

xE 1 \yX 1 (x 1 ),x l ]}, (1.42) 

where the universal number A/ was found to be A/ = 1 — 0.188/AT + 0(1/N 2 ). 

In Section |V| we will present the details of our computations of the reduced scaling 
functions of the renormalized-classical regime in a 1/N expansion. Our results agree with 
those of Chakravarty et. a/0; however we are also able to obtain a number of universal 
amplitudes which were previously only determined by numerical simulations. 

Another of our new results here is the low T dependence {x% <C 1) of the scaling functions 
for the uniform susceptibility and specific heat: 

fii(xi) 



tixi Nit 

V 1 (x 1 ) = N-l (1.43) 

These results have also been obtained independently by Hasenfratz et. a/.@ The first of these 
results will be important to us later in the comparison with experiments. The second result 
simply implies that the number of low energy degrees of freedom are the N — 1 spin waves. 
The two results together also imply that the Wilson ration W ~ 1/T at low temperatures 
provided g < g c . 



G. Quantum disordered region (x2 -C 1) 

At T = (x 2 = 0), the ground state of the antiferromagnet has a gap towards all 
excitations in this region. Unlike both previous regions, therefore, finite T is almost always 
a weak perturbation on the T = results; all finite temperature corrections are accompanied 
by factors of exp(— A/(&bT)) = exp(— l/a^) <C 1. Furthermore, the ground state is rather 
well described by the N = oo theory. We will therefore refrain, here, from giving complete 
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expressions for all the observables; we refer the reader to Section |III A 2| for the exact results 
at iV = oo. 

However, thermal effects and 1/N corrections are important at measurement frequencies 
smaller than A; in this region the dynamics is controlled by a dilute concentration of ther- 
mally excited quasiparticles. The dissipative effects of such quasiparticles will be discussed 
in Section |[V|. 

We now discuss the plan of the remainder of the paper. In Section |TJ we present a 
phenomenological derivation of the scaling forms used above. Section [III] introduces the 
quantum O(N) non-linear sigma model and presents the complete solution for all the scaling 
functions at N = oo. The formal structure of the model at order 1/N is also discussed. The 
remainder of the paper contains the details of the calculations. The calculations are discussed 
for the quantum-critical, renormalized-classical, and quantum-disordered regions in turn in 
Sections [IV], [V], [VTj. Finally in Section [VI I| the comparison with experimental results is 
presented. The appendixes contain a discussion of the effects of disorder and Berry phases, 
results of a Monte-Carlo simulation, and some technical details. 



II. PHENOMENOLOGICAL DERIVATION OF SCALING FORMS 



In this section we will present a phenomenological derivation of the scaling forms ( |1.10a ) 



( |1.19a| ) for the order-parameter dynamic susceptibility, the scaling forms (|1.10b|) , (|1.19b ) for 



the uniform spin susceptibility, and the scaling forms ( |1.10c| ), ( |1.19cj ) for the specific heat. 
These are valid in the vicinity of a quantum phase transition in a two-dimensional quantum 
Heisenberg antiferromagnet from a state with long-range Neel order to a spin-fluid state. We 
will only consider the case in which the Neel order parameter is a 3-vector. The following 
discussion does not explicitly refer to the quantum non-linear sigma model. It should instead 
be regarded more generally as a study of the consequences of the scaling hypothesis on a 
quantum phase transition in a Heisenberg antiferromagnet. The non-linear sigma model 
provides a realization and verification of these hypothesis for a particular field theory. 
Let us first present the precise ingredients from which our results follow. 

1. The spin-wave velocity, c, should be non-singular at the T = 0, quantum fixed-point 
separating the two phases. We will also assume, for simplicity, that there is no spatial 
anisotropy; this assumption is not crucial and our results can be easily extended to 
include quantum transitions in anisotropic systems like spin chains coupled in a plane. 
These systems will of course have two spin-wave velocities, whose effects can be ab- 
sorbed into a rescaling of lengths. In the presence of spatial isotropy, and for the case 
of the vector order-parameter, the non- singularity of the spin-wave velocity implies 
that the critical field-theory has the Lorentz invariance of 2+1 dimensions. 

2. The antiferromagnet in the vicinity of the critical point satisfies 'hyperscaling'S hy- 
pothesis. For a quantum transition in D = 2 + 1 dimensions, this hypothesis implies 
that the singular part of the free energy density (J-'g) at T = has the form 

T s = hcTC D (2.1) 
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where £ is the correlation length which diverges at the transition (of course, on the 
ordered side £ = The number T is dimensionless, and like all such numbers at 
the critical point, it is expected to be universal. The statement of the universality of 
T is known in the literature as the hypothesis of 'two-scale factor' universalityii'B. 

3. Turning on a finite temperature places the critical field theory in a slab geometry 
which is infinite in the two spatial directions, but of finite length, 

Lr = £=, (2.2) 
k B T 

in the imaginary time (r) direction. The consequences of a finite T can therefore be 
deduced by the principles of finite-size scalingil. 

In the following we measure temperature, time, and length in units which make ks, fo, 
and c equal to unity. 



A. Staggered susceptibility, g < g c 

Consider first the application of the scaling hypothesis to the staggered spin susceptibility 
for g < g c . A straightforward application of finite-size scaling yields 

X.(h,u) = AL 2 -^ (kL T ,uL T , |^ 

4 -$(k,uJ,^-), (2.3) 



T 2 -v V Lr 

where A is a non-universal amplitude, and the scaling function <3> is universal upto an 
overall pre-factor. In particular there are no non-universal metric factors^ in any of the 
three arguments of $. We now wish to eliminate the dependence of this result on £j. It has 
been argued recentlyHJll that the result ( |2.1| ) can be extended to deduce a simple, universal 
relationship between p s and £j, valid in the limit £j — > oo: 



hcT 
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(2.4) 



where T is another universal number. Now, from (|1.13|) , ( |2.4j) and (2^2) we see that that 
the third argument 

f = ^ = ^, (2.5) 

L T Ps 27T 

We therefore define a new universal function $i s 

~ /_ NT \ 

$i a (k,u 1 x 1 ) = $\k,u 1 —x 1 \. (2.6) 
We shall soon see that $i s is the same function as in ( |1.10a| ). 
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To make $i s completely universal, we have to fix its overall scale, which we now do. 
First, we notice that in the renormalized classical region (x\ <C 1), there is a Goldstone 
regime (k£j < 1, fc( > 1, see Fig. |2] and Section [TFD of non-interacting spin-waves. In 
this regime, the hydrodynamics predicts that the static staggered susceptibility has a simple 
form: 



N 2 



NJ p s k 2 ' 

x x < 1, k£j < 1, H > 1- (2.7) 

The scale of the susceptibility has been set by the magnitude of T = staggered moment. 
The factor of (1 — 1/N) arises from the fact that only N —1 transverse modes contribute the 
Goldstone singularity, while the longitudinal mode is massive; after rotationally averaging 
this induces a factor of (N—l)/N. We now demand that in the appropriate limit, the scaling- 



form ( |2.3| ) obey this Goldstone form. The key constraint is that (|2.7|) is T independent. It 
is easy to show that this can be satisfied by (|2.3| ) only if 

- _ b 

$ u (k,u = 0,xi)= =2— ; 

k x{ 

Xl < 1, kx x < 1, ke N/ ^ N - 2)xi) > 1, (2.8) 

for some constant b. In the last restriction on k we have used ( |1.40|) , valid for the O(N) 
sigma model and neglected pre-exponential factors of x\. The prefactor in ( |2.8j ) is of course 
arbitrary. We now make the specific choice, b = 1 — 1/N and thus specify the overall scale 
of $ ls . Comparing Q2.3|), (^.6[), (^.7[), and (|2.8|), we can now fix the value of the prefactor A 

A^(JL)\ (2.9) 
p s \2icp s J 

Inserting this value of A into ( |2.3| ), we obtain ( |1.10a| ) as desired. 



B. Staggered susceptibility, g > g c 

A similar analysis can be carried out in the spin-fluid phase for the staggered suscepti- 



bility and its scaling function $2s as in (1.19a). We remind the reader that our theory is 



valid only for antiferromagnets with a 3-vector order parameter, in which case the quan- 
tum disordered phase is expected to have only integer spin excitationsiBEJ; in particular if 
spin- 1/2 spinons are present at intermediate scales, they are always confined at the longest 
distances. 

At T = the equal-time order parameter correlation function will decay in space with 
a correlation length £. As the theory has a Lorentz invariance, this implies that the gap 
towards spin-1 excitations, A, is 

he 

A = -, (2.10) 
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where we have momentarily reinserted explicit factors of % and c. An application of finite- 
size scaling, very similar to that for the ordered side, now yields immediately the scaling 
function (we now return to units in which k% = h = c = 1) 

Xs(k,u) = 7 ^$ 2s (k,u!,^). (2.11) 



The only non-universal components on the right-hand-side are the amplitude A and the 
related overall scale of the function § 2s . As before, we will fix this scale by matching 
with an experimental observable at T = 0. Let us therefore think about the nature of 
the spectrum in the spin-fluid phase at T = 0. As all excitations have a gap, the spin 1 
quasiparticle should have an infinite lifetime for energies close enough to the threshold A 
(this quasiparticle appears as the bound state of two spinons in large M theories of SU (M) 
and Sp(M) antiferromagnetsEllli). Further, the Lorentz invariance of the theory implies 
that the dispersion spectrum, uj^ of this spin-1 quasiparticle is given by 



uj k = Vk 2 + A 2 (2.12) 

for small enough k. These facts combine to imply the following form for Xs at T = 

A 



k 2 -(u + ief + A 2 ' 

T = ,A;<A, |tu-A|<A, (2.13) 



where e is a positive infinitesimal and A is the spin-1 quasiparticle residue. This residue 
can be experimentally measured by examining the imaginary part of the local susceptibility, 
\mxL- Using ( |1.16| ) and ( |2.13| ) we find 

A 

IhixlM = —Q{uj - A) ; T = 0, \uj - A| < A. (2.14) 

The discontinuity in the local dynamic susceptibility is therefore precisely a quarter of the 
spin-1 quasiparticle amplitude. We now demand that the scaling form ( 2.11 ) satisfy ( |2.13| ) 
as T — > 0. A little experimentation shows that this is only possible if 

2-T) 

§ 2s (k,u},x 2 ) 



(kx 2 ) 2 - (ujx 2 + is) 2 + 1 
x 2 < 1, kx 2 < 1, \px 2 - 1| < 1. (2.15) 

Here we have arbitrarily set the overall scale of Q 2s ; this function is now completely universal. 
Finally, the derivation of ( 1.19a ) is completed by obtaining the amplitude A in ( f2.11| ): 

* = < 2 - 16 > 

The prefactor A is expected to be non-singular as g \ g c ; therefore the quasiparticle ampli- 
tude A must vanish as 

A~(g-g c )o". (2.17) 

This was noted earlier in ( |1.18| ). 
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C. Uniform susceptibility 



The key ingredient in the determination of the scaling form for Xu is the realization that 
Xu is simply a stiffness related to twists in boundary conditions on the system along the 
imaginary time (r) direction^. A uniform magnetic field on the antiferromagnet causes a 
precession of all the spins at the same rate. The relative angle between any two spins remains 
unchanged. Therefore by transforming to a rotating reference frame almost all vestiges of 
the magnetic field can be removed. However the partition function in the laboratory frame 
had periodic boundary conditions along the r direction, implying that the system in the 
rotating frame has a twist in its boundary condition. The susceptibility is the response to 
such a twist, which is precisely the spin stiffness. 

The Lorentz invariance of the theory now implies that the scaling properties of Xu should 
be the same as those of p s , which is the stiffness for twists about the spatial boundary 
conditions. In other words, the scaling dimension of Xu is exactly D — 2, and at T = 
the combination approaches a universal number as g /* g c . The scaling laws (|1.10b ), 
( |1.19b|) are now a completely straightforward consequence of the principles of finite-size 
scaling. Unlike Xs, there is no need for any normalization condition to set the overall scale 
of the scaling function^ It is already fixed to a universal value by the hypothesis of two-scale 
factor universality 5,5 



D. Specific heat 

Consider the D = 2 + 1 dimensional Lorentz invariant theory in a slab geometry in the 
vicinity of g = g c . An early paper by Fisher and de Gennesll argued by extending ( |2.1| ) to 
finite sizes, that the free energy density T must have the following dependence on the size, 
L T of the finite dimension 



^ = ^o + T7^> \¥- > (2-18) 




where JF is the bulk free energy density and we have assumed that g < g c . The function 
tp is universal at all x; there are no non-universal metric factors in either the argument or 
the scale of yS. The scaling function ( p.. 10c ) for Cy now follows immediately from the 



thermodynamic relationship between T and Cy, the relationship ( j2.4|) between p s and £j, 
and the relationship ( p.2|) between L T and T. An entirely analogous argument can be made 
for g > g c . 

We note that a related result has been discussed recently by Castro Neto and Fradkina 
in the context of 2+1 dimensional quantum systems; they have also discussed an interesting 
connection to the Zamalodchikov's C-theorenS We have chosen a numerical prefactor of 
the scaling function ^ in ( |1.10c| ) which is the specific heat of a single gapless bose degree of 
freedom with dispersion u = ck in 2 dimensions. The number ^i(T — ► 0) is thus a measure 
of the effective number of such modes in the ground state. For g < g c , this number should 
equal N — 1, the number of spin- wave modes in the ordered state. For g = g c , the number 
\I/i(oo) is probably irrational and will be calculated later in this paper to order 1/N. Finally, 
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in the quantum disordered phase, g > g c , there are no gapless modes and we should have 
* 2 (0) = 0. 

A very similar connection between the effective number of gapless modes and Cy was 
established some time ago for 1 + 1 dimensional spin chainslll; in this case Cy ~ T with 
Cy I T universal and related to the central charge of a conformal field theory, which in effect 
measures the number of gapless quasiparticle modes. 



III. THE QUANTUM O(N) NON-LINEAR SIGMA MODEL 

In this section we will discuss the O(N) quantum non-linear sigma model field theory. 
This theory may be viewed as the simplest model which displays a quantum phase transition 
in 2+1 dimensions. Moreover, a microscopic connection between weakly frustrated antifer- 
romagnets with short-range Neel order and the 0(3) sigma model can also be established!. 

There are several subtle and difficult questions relating to the consequences of Berry 
phase terms which are present in the antiferromagnet but are absent in the sigma modellll. In 
this paper we will simply neglect the effects of the Berry phase terms. In Appendix |B] we will 
present some circumstantial evidence, in computations for SU (M) antiferromagnets, which 
suggests that these Berry phases are irrelevant both in the Neel phase and at the quantum 
critical point, but do significantly modify the properties of the quantum disordered phase. 
In critical phenomena terminology this implies that the Berry phases are 'dangerously- 
irrelevant'. As all of our scaling functions are properties of flows in the vicinity of the 
quantum-critical point, we do not expect any modifications of our results by Berry phase 
effects in this scenario. Berry phases will however modify the corrections to scaling. 

In passing, we note that there is an alternative expansion which could have been used to 
obtain the scaling functions of this paper: this is the large M expansion about antiferromag- 
nets with SU (M) or Sp(M) symmetry0@^i. However, the presence of a gapless gauge field 
in the perturbative 1/M corrections, makes this expansion somewhat more involved than 
the 0(N) expansion. For general M, N, both approaches predict that the lowest-lying non- 
zero spin excitation above the quantum-disordered ground state in an antiferromagnet with 
collinear orderly carries spin S = 1; the detailed structure of the spectrum at higher energies 
is however different in the two theories. The results of this paper show that the 0(N) expan- 
sion is numerically much more accurate at N = 3 than is the Sp(M) or SU (M) expansion 
at M = 2. This can be seen immediately by comparing the values of rj in the two theories: 
the large M theory gives r\ — 1 — 0(1/ M) while the large N theory has rj = + 0(1/ N) 
- compare this with the known valued in the D = 3 classical Heisenberg model t] « 0.028. 
We shall see below that, at N = 3, the 1/N corrections in the O(N) model to the universal 
scaling functions are almost always less than about 20% of the leading N = oo term. 

We will begin this section by a definition of the quantum O(N) non- linear sigma model. 
The first subsection will present its exact solution at N = oo. This solution has also been 
discussed earlier by Rosenstein et. alM, although they did not emphasize the universal 
scaling properties of the solution. For clarity, we will repeat some of the step in Refill 
and will then explicitly compute all of the scaling functions introduced in Section [I] in the 
N = oo limit. The N = oo computations are also similar to earlier studies of finite-size 
scaling properties of the spherical modeS - in our case the inverse temperature plays the 
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role of the finite-size along the time direction. The second subsection will present complete 
formal expressions for the staggered and uniform susceptibility which are correct to order 
1/N; to the best of our knowledge, these constitute the first computations of finite-size 
corrections to two-loop order in any system. Subsequent sections will manipulate these 
expressions into the appropriate scaling forms for the three regions of Fig. |I[ The structure 
of the 1/N corrections is however rather involved and the casual reader may be satisfied by 
studying only the N = oo solution of Section [111 A|. Even this limited solution is quite rich 



and instructive; its main shortcoming is the absence of spin-wave damping and anomalous 
dimensions which appear only at order 1/N. 

The O(N) non-linear sigma model is defined by the functional integral 



Z = J Vn E 5(nj - 1) exp ( J d 2 r j <h 



(V r n £ ) 2 + \(d T n £ ) 
Co 



(3.1) 



where the index i runs from 1 to N, p° s is the bare spin stiffness, and Cq is the bare spin-wave 
velocity. Both p° and Cq differ from their renormalized values p s and c; however at iV = oo 
we will find Cq = c, although the renormalization of Cq will be quite crucial in subsequent 
sections. Further analysis is simply expressed in terms of the coupling constant, 
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Nhc 



(3.2) 



which has the units of inverse-length. We will find that the quantum transition occurs at a 
g of order unity. This implies that we have to choose p° ~ N in the large N limit. In the 
remainder of this section we will use units such that ks = h = Cq = 1. The large N analysis 
of Z begins with the introduction of the rescaled field, 



Nn e 



(3.3) 



and the imposition of the constraint by a Lagrange multiplier A. This transforms Z into 



1 



d 2 r / dr 
Jo 



(y T ntY + (d T n t y + i\(hf - N) 



(3.4) 



This action is quadratic in the ne, which can therefore be integrated out. This induces an 
effective action for the A field which has the useful feature of having all its iV dependence in a 
prefactor. Therefore, for large N the A functional integral can be evaluated by using its value 
at its saddle-point. Terms higher-order in 1/N can be obtained by a systematic expansion 
of the functional integral about this saddle-point. We parametrize the saddle-point value of 
A by 



m 



(3.5) 



where the 'mass' m is to be determined by solving the constraint equation n\ = 1, order- 
by-order in 1/N. This value of m can then be used to obtain a 1/N expansion of the n^-n^ 
correlator and hence of all the observables related to the staggered susceptibility. 

Determination of the uniform susceptibility requires introduction of a slowly-varying 
magnetic field B into the non-linear sigma model. In the 0(3) model such a field causes 
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a precession of the local order parameter about the magnetic field axis. This precession is 
realized in the 0(3) sigma model by the substitution!! 

d T n — > d T n — i^^-B x n (3.6) 

in functional integral Z\ here g^s/his the Bohr magneton for a single spin; in the remainder 
of this section we will measure the field in units of gfis/k and hence omit explicit factors of 
gUs/h. For the general 0(N), the analog of the magnetic field is a second-rank antisymmet- 
ric tensor bi m which causes a precession of the spin- components lying in the plane defined 
by the £, m directions. For iV = 3, bg m is related to B by 

bim — timpBp. (3-7) 

The full action in the presence of the b field is defined by ( |3.1|) and the substitution 

d T rii -> d T n £ - ib em n m . (3.8) 

The uniform susceptibility, Xu is now obtained by evaluating In Z in powers of b and picking 
out the coefficient of the quadratic term: 



1 Cp 1 

X "=2T^C l0gZ ' (3 - 9) 



where V is the volume of the system. 



A. Solution at N = oo 

The 0(N) sigma model can be solved exactly at N = oo. Closed form expressions for 
all the scaling functions introduced in Section | can be easily obtained. We begin with the 
staggered spin susceptibility which is given by 

Xs{k,u)) = ?^-— — J_ (3.10) 

We have introduced here the mass m which is the value of m at iV = oo, and S = Z$S which 
is a rescaling factor between the actual susceptibility of a quantum spin-S antiferromagnet 
and the susceptibility of the unit n-field in the 0(3) sigma model. The renormalization 
factor Zs accounts for the fluctuations at short scales, which have to be integrated out in 
the derivation of the sigma model from the original spin Hamiltonian. Expressions for all 
the experimental observables dependent upon Xs can now be easily obtained; we will refrain 
from giving explicit expressions. One quantity we will need is the correlation length £, which 
we defined earlier from the long-distance decay ~ e~ r ^ of the equal-time n-n correlation 
function. Such a decay is present at all finite T for all values of g. At iV = oo we find 




(3.11) 
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Now turn to the uniform spin susceptibility. There is no damping at N = oo and hence 
the spin-diffusion constant is infinite. We therefore consider only the value of the static 
susceptibility. Evaluating (|3.9|) at iV = oo we find 



where = k 2 + tuq, and u n is the Matsubara frequency which takes the values 2irnT with 
n integer. The frequency summation and the subsequent momentum integration can be 
performed exactly. Moreover, the momentum integration is convergent in the ultraviolet, 
and final result depends only on rriQ and T: 



T 

y"t = _ 
7T 



m /T 

H^— wfe mo/T -r 



(3.13) 



Note the absence of direct dependence on the coupling g. 

Consider next the free energy density, J 7 , from which Cy can be obtained by taking two 
temperature derivatives. At iV = oo iSH Fcan be directly obtained from the saddle-point value 
of the effective action. The result is 



NT ^ r d 2 k 2 2 2n Nm2 
-^=^-Ey t-^ Mfc 2 + ^ + mg) 

= ^/0Ml-e~^) 
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A r d 2 kdu , /l 9 9 9 . Am^ 

+ t/^ 10 ^ +w +m »'-^ (3 - 14 » 

It now remains to determine the dependence of mo on g and T. The constraint n 2 a = 1 
takes the following form at M = oo: 

T ^J 4^ki + S n + m 2 =L (3 - 15) 

It is easy to see that the momentum summation is divergent in the ultraviolet, and it is 
therefore necessary to introduce an cut-off. We use a relativistic Pauli-Villars cut-off at the 
momentum scale A, assumed to be much larger than the temperature. This transforms the 
constraint equation into 

roo d 2 k ( g g \ 

-Wo 4^ \k 2 + u 2 + m 2 ~ k 2 + u 2 + A 2 J = ' (3 ' 16) 

The momentum integration is now convergent in the ultraviolet and can be performed ex- 
actly. The subsequent frequency summation is also tractable and yields 

9 _T / s inh(A/2rn = 

2tt \sinh(m /2T) J v ; 

Finally, we can solve for the dependence of m on T, g, and A 



24 



m 



2Tarcsinh exp 



| sinh ( — 

gTJ V2T, 



(3.18) 



In the limit T < A, this equation can be rewritten as 



mo = 2Tarcsinh 



1 / 2tt 1 V 



where 



9c 



47T 



(3.19) 



(3.20) 



By examining the T — > limit of this equation, it is immediately apparent that the behavior 
of £ is quite different depending upon whether g is smaller, larger, or close to the critical 
value g c . For g > g c , £ approaches a finite value as T — > 0, while for g < g c , £ diverges as 
T — > 0. The finite, T = asymptote of £ for g > g c itself diverges as g approaches g c . We 
find 



Z(T = 0,g>g c )~(g-g c )- 1 . 
This identifies the N = oo value of the exponent v to be 

i/ = 1. 



(3.21) 



(3.22) 



For g < g c these is a Josephson scaleo £j which diverges with the same exponent v as g 
approaches g c . However, this scale is not present in the N = oo theory: this is because the 
exponent rj = at N = oo, making the critical and Goldstone fluctuations indistinguishable. 
At g = g c , we have from 



where the number 



9 = 2 log 



m = 6T, 



'v^+l^ 



0.962424 



(3.23) 



(3.24) 



will occur frequently in our analysis. We see therefore that the correlation length scales with 
1/T at the critical point. 

We now examine the scaling limit of the above results. This limit is obtained when the 
temperature T, and the deviation from criticality \g — g c \/g1 are both much smaller than 
the upper cutoff A. We will consider the cases of g smaller or greater than g c separately. 



1. Scaling properties for g < g c 

From the discussions in Sections |IA| and it is clear that the results become simple after 
they have been expressed in terms of the T = ordered staggered moment N and the fully 



25 



renormalized, T = spin-stiffness p s . In Appendix |D] we have performed a 1/N expansion 
of these T = variables. At iV = oo their exact values are 



Ps 


= "( 


\9 


0c/ 


N 2 


= S 2 


(- 


-) 






9c J 



(3.25) 



We now use these results with the expression ( |3.18| ) for £ to study the limit p s /iV,T < A. It 
is not difficult to show then that £ takes the scaling form ( 1.30| ) where the scaling function 
X x is 

Xl(zi) = 2 arcsinh f — — J . (3.26) 

Thus, in the renormalized classical region X\ <C 1, we have 

X x = e- 1/a;i ; Xl < 1, (3.27) 

implying a correlation length which is exponentially large, while in the quantum-critical 
region (an ^> 1) we find 

X, = 2 log - ^ ; an » i. ( 3 - 28 ) 

indicating that the correlation length, £ = l/(TXi) scales with 1/T. 

The scaling functions for all the other observables now follow in a straightforward manner. 
The scaling function <3> ls for the staggered susceptibility in ( |1.10a| ) and the exponent rj can 
be deduced from Q), (^TJ, and ([DSD} to be 

77 = 0, 

57, *i) = =3— -L (3.29) 

fc - (cJ + ^) 2 + X 1 2 (xi) 

The scaling functions for the observables dependent upon Xs now follow immediately. For 
the structure factor S(k) we have the scaling function Hi in ( [1.281 ) which is 

E 1 (k,x 1 ) = — 2 coth ^ 2 + Xl2(Xl)]1/2 . (3.30) 

U ' 2[k 2 + Xl{ Xl )Y^ 2 1 ; 

The local susceptibility Imxi(^) obeys (|1.31|) with the scaling function F\ given by 

^(57, x x ) = i [9(cJ - Xx(xi)) - 9(-U - Xx(xi))] , (3.31) 

where 6{x) is the unit step function. The presence of a gap in F\ at finite T is an artifact 
of the N = oo theory and will be cured upon including 1/N corrections. A consequence of 
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this gap is that there is no relaxation of nuclear spins at N = oo and the scaling function 
for 1/Ti in (|L35j) is 



R^xi) = 0. 



(3.32) 



The result for the uniform susceptibility, xti ( |3. 13|) can also be collapsed into the scaling 
form ( |1.24j ). We use ( |3.13| ), the expression ( |3.18j ) for m , and ( |3.25| ) to obtain 



1 



V4 + e- 2 /^i , fe~ 1/xi 
-arcsinh 



ttxi 7re l / Xl \ 2 J 

In the renormalized classical limit x\ <C 1 this function has the limiting behavior 



fii(xi) 



1 1 
— + - 



TTXi IT 

while in the quantum-critical region (x\ 3> 1) it obeys 



Xi <C 1, 



fii(xi) 



7T 



log 



m 


i 4 ' 




5xi. 



Xi > 1. 



(3.33) 



(3.34) 



(3.35) 



We turn now to the free-energy density, T and the specific heat Cy. We will evaluate 
J-(T) —TiT = 0) at a fixed value of g. The calculations are performed most easily if we use 
relativistic cutoff and write the value of g c in the following form 



d 3 P 1 

8^3 P2 ; 



(3.36) 



where P = (k, u) is the relativistic 3-momentum and the integral is suitably regulated in 
the ultraviolet. We now apply fl3.25|) and the result mo(T = 0) = to fl3.14j) to obtain 



^(T)-^(0) = NTj^ log(l - e~^ T ) 



+- 



N rd 3 P 



in" 



log 



P 2 + ml 

P 2 







m 2 ' 




mlp s 


P 2 




2 



(3.37) 



Using m = TXi(xi), and simplifying the above integrals, we get 
JF(T) - ^(0) = NT 3 



oo 

— log(l — e" 

Xi(a:i) 27T 



Xf^i) X 2 { Xl ) 



12tt 



Attxi 



(3.38) 



Finally we use the thermodynamic relation Cy = —Td 2 TjdT 2 to obtain the scaling function 
\I / i(xi) for the specific heat as defined in ( |1 . 10c| ) 



A% d 2 
3^(3)xi (ixj 



.r, 



— log 1 - e 

Xi(xi) 2tt 



127T 4lTXi 
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(3.39) 



We recall that the function Xi(xi) is given in ( |3.26| ). It is clearly that the universal crossover 
function is quite non-trivial even at iV = oo. In the renormalized classical limit, 

X\ — > it has the value 



*i(0) = N, 



(3.40) 



which is the number of gapless spin wave modes (upto relative order 1/N). The normaliza- 
tion in ( |1.10c|) was chosen to make this result simple. In the quantum-critical limit, X\ — > oo, 
the integrals cannot be analytically evaluated. However, one of us has recently shownEil that 
application of some unusual identities of polylogarithmic functions can be used to prove that 
the integrals reduce a surprisingly simple result 



AN 

*i(oo) = — . 



(3.41) 



We have no physical understanding of why this number is rational. 

The reduced scaling functions of the renormalized classical region, describing the 
crossover from Goldstone to classical, thermal disorder (Fig. [|) can also be easily obtained. 
From ( |1.41| ), ( |1.42| ), and ( |3.30| ) we find the scaling function f(y) for the structure factor 
S(k) 



f(y) 



y 2 + i 



(3.42) 



2. Scaling properties for g > g c 

The first step is to obtain the value of the T = gap. From (|3.19|) we have the exact 
N = oo result 

A = 4tt ( - - - ] . (3.43) 
\9c 9 J 

Second, we also need the spin-1 quasiparticle amplitude A. From ( |3.10| ) we obtain immedi- 
ately at N = oo 

A= 9 -^. (3.44) 

From ( 3.18|) we can now deduce the scaling function X2 as defined in ( |1.30 ) for the 
correlation length 

/ l/(2:z 2 ) \ 

X 2 (x 2 ) = 2 arcsinh f — - — I . (3.45) 

(Recall that x 2 = T/A.) This function has the following asymptotic limits in the quantum 
disordered and quantum-critical regions respectively 
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X 2 {x 2 ) 



— + 2e~ 1/x2 ; 
x 2 



x 2 < 1, 
x 2 > 1. 



(3.46) 



These results imply a correlation length which is of order A 1 and T 1 respectively. 

The scaling function fl 2 (x 2 ) f° r the uniform static susceptibility in ( 1.24 ) can be obtained 
from (|3~T3|) combined with (gig) and Q). We find 



tt 2 (x 2 ) 



1 V4 + e 1 /^ 

H 777^- s ar csinh 



,1/(2*2)' 



a/(2x 2 ) 



27nr 2 Tre^ 

In the quantum-disordered region (x 2 <C 1), this function is exponentially small 

g-l/32 



^2(^2) 



vrx 2 



x 2 < 1, 



while in the quantum-critical region it implies a Xu °f order the temperature 



^2(^2) = ^7 log 1 



1 



2 

5x 2 



x 2 > 1. 



(3.47) 



(3.48) 



(3.49) 



Next we consider the free-energy density, T and the specific heat Cy. The following 
relation between the T = gap, A, and the coupling g will be useful 

d 3 P 



1 

9 



1 



(3.50) 



^3 p2 + A 2 ' 

where P = (k, to) is the relativistic 3-momentum. Note that we again use relativistic cutoff 
for the value of g c . Applying this result to ( |3.14| ), we obtain 

d 2 k 



HT)~m = NT /^log(l 



+ - 



N rd 3 P 



log 



'P 2 + ml\ mg-A 5 



P 2 + A 2 / P 2 + A 2 



(3.51) 



Using m = TX 2 (x 2 ), and simplifying the above integrals, we get 
JF(T) - ^(0) = NT 3 



77-l0g(l 

X 2 (a; 2 ) 27T 



Xf(x 2 ) , 3x 2 Xf(x 2 ) - 1 



12tt 



24ttx3 



(3.52) 



Finally we use the thermodynamic relation Cy = —Td 2 TjdT 2 to obtain the scaling function 
^2(^2) for the specific heat as defined in ( |1.19c|) 



^2{x 2 ) 



Nn d 2 
3C(3)x 2 dx\ 



.1, ■-. 



— log(l - e 

X 2 (x 2 ) 2,71 



X 3 (x 2 ) 3x 2 X 2 (x 2 ) - 1 



12tt 



2A7TX 3 , 



(3.53) 
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where the function X^ix-i) is given in ([3.45 ) . In the quantum disordered limit, x 2 — > 0, we 
have ^2(0) = corresponding to the absence of any gapless degrees of freedom. Recall also 
that ^2(00) = *i(oo). 

The scaling functions for the remaining observables are identical to those obtained above 
for g < g c after the substitution Xi(xx) — > X 2 (x2). Thus e.g. the scaling function $ 2s for 
Xs in (|1.19a ) differs from the function <3>i s in ( p.29|) only by this substitution. 



B. 1/N corrections 

We now present formal expressions for the modifications to the uniform and staggered 
spin susceptibilities at order 1/N. These expressions will be used in the following sec- 
tions to determine the scaling properties of the quantum-critical, renormalized-classical and 
quantum-disordered regions. 

The formal structure of the 1/N expansion has been reviewed in the book by Polyakovil 
(see also E1§1). The leading corrections arise from considering the fluctuations of the A field 
about its saddle point. The contribution to the staggered susceptibility, or equivalently the 
propagator of the field, is given by the Feynman graph in Fig. |] This leads immediately 
to 

gS 2 1 

Xs(k, iu n ) = ^ fc2 + ^ + m2 + s(M ^ r ( 3 - 54 ) 

where £ is the self-energy arising from the A fluctuations. It is convenient to absorb the 
value of £(0,0) into the mass m 2 ; all our expressions for £ will therefore always include a 
subtraction which makes £(0,0) = 0. The leading contribution to £ is 

E(k,iuj n ) = 

2 r d 2 g G (k + q, iu n + ie n ) - G (q, ie n ) 

where 1 /II is the propagator of the A field 
U(q,ie n ) = 

^" G o(q + q.i,ie n + iCLn)G (q 1 ,iCl n ), (3.56) 
and Go is the propagator of the ne field at iV = 00 

Go{k '™» ) = k* + (3 ' 57) 

It now remains to determine m 2 . The value of m 2 is set by solving the constraint equation 
n 2 = 1, or 

r E/0X.(*,*O = ^, (3-58) 
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order by order in 1/N. At iV = oo, the dependence of mo on g and T is determined by 
solving the following equation exactly 



(3.59) 



This has of course already been done in the previous subsection. We represent the 1/N 
corrections to m 2 by 5m 2 , where m 2 = m\ + 5m 2 . Upon examining the 1/N corrections to 
( j3.58| ) we find in a straightforward manner 



5m 



Ki{T, mo) 



K 2 (T, m )' 
where the constants K\, K 2 are 

f d 2 k 

fTi(T,m )=T^ / — E(k,iu n )Gl(k,iuj n ) 
K 2 (T,m )=TY: J 0G 2 (Mc n ) 



^3.60) 



87rm 



coth 



mo 
2T 



(3.61) 



Note that K\, K 2 depend only upon T and m , and do not depend directly upon the value 
of the coupling g. Comparing ( |3.61| ) with (|3.55|) , (|3.56| ) and ( |3.57|) , the expression for K\ 
can be manipulated into the following 



Ki(T,m ) = ~T y £,J 
/ 1 dll(q,ie 



d 2 q 



4mo dmo 



Air 2 U(q,ie n ) 
+ G {q,ie n )K 2 (m ,T) 



(3.62) 



These results for m 2 , K\ and K 2 will be quite useful in subsequent sections. 

Expressions for observables that depend upon \s can now be obtained as before. In 
particular, the equal-time staggered structure factor S(k) is given by 



s (k) = Tj2Xs(k,iuj n ). 

Using ( 3.54j ) we see that to order 1/N this can be rewritten as 
S(k) _ [ E^(£m 2 + Z(k,tu n ))G 2 (k,iu n ) 



(3.63) 



S (k) 



S (k) 



(3.64) 



where S (k) is the structure factor at iV = oo 

9~S 2 



S (k) 



N 



(3.65) 
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The correlation length, £, is defined by the pole of S(k) in the complex k plane which is 
closest to the real k axis; this pole is at k = From (|3.63| ) it is clear that, at any finite 
T, the poles of S(k) are simply the poles of Xs(k,iuj n ) for all u n . It is also clear that the 
pole closest to the real-axis is that associated with Xs{k-,iu n = 0). The location of this pole 
can be simply determined from ( |3.54j ) in a 1/N expansion. We find 



£ 1 = m + - — E(/c = im , iu n = 0). (3.66) 

ZTTIq 

Finally, the residue of the pole in the structure factor, which sets the overall scale for 
exponentially decaying correlations, can be determined in a similar manner: 

dT,(k = imo,iu n = 0)\ 1 




dk 2 k 2 + i- 



k close to i£~ 



(3.67) 



The results for £ and S(k) will be of great use to us in the subsequent sections. 

Now we turn to a consideration of the 1/N corrections to the static uniform susceptibility, 
Xu- This is a four-point function of the field and the Feynman graphs which contribute 
at order 1/N are shown in Fig. The evaluation of these graphs is, in principle, quite 
straightforward, but rather tedious. After some fairly lengthy frequency summations, we 
obtained the following, surprisingly compact result: 



d 2 k 

4^ 



777 



-2n',-^(m 2 -m2) 



- K 3 (T, mo), 



(3.68) 



where the energy = (k 2 



m o) 1//2 was introduced earlier, n k = ?7(efc) is the Bose function 



me 



(3.69) 



and n' k = dn(ek) / de^, etc. (note that the symbol n is used both for the Bose function and 
the non-linear sigma model field - the appropriate meaning should however be clear from 
the context). The very first term in ( |3.68| ) is identical to the N = oo result of ( |3.12| ), while 
the second term arises from effective mass renormalization of the the N = oo graphs. The 
remaining 1/N corrections are contained in the function K 3 (T,m ) which is given by 



K 3 (T, mo) 



-TY 
N ^ 



~'k+q 



d 2 q 



Air 2 U(q, iu n ) J An 2 e k 



d 2 kn'l 



+ 4 



+ uo 2 ) 2 -Ae 2 e 2 



k c k+q 



q 1 n 



(3.70) 



We will evaluate this expression in the later sections on the different critical regimes. 

Lastly, the 1/N corrections to the free energy density, T. The 1/N corrections to ( 3.14|) 
arise from the functional determinant of the integration over the A field. The propagator 
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of the A field is II as defined in ( |3.56| ) and the correction to the free energy in simply 
(l/2)Trlog!l. We have therefore 




d 2 k r 
4^ . 



\og{k 2 + u 2 n + m 2 ) 



+ T7log(n(Mu; n )) 



Nml 



(3.71) 



■2fJ 



Note that we have not included the 1/N correction to the mass i(X) = m. This is because 
by construction dJ r /d{\) = at iV = oo. This immediately implies that the correction to 
(A) at order 1/N will modify T only at relative order 1/N 2 . Thus may as well use i(X) = m 
to compute T to order 1/N . 



This section will present expressions for the scaling functions in the quantum critical 
region to order 1/N. We will restrict our discussion here to the critical point at g = g c . 
This point will be accessed by taking the X\ — > oo limit from the ordered side. We will thus 
present explicit results for the scaling functions $i s (A;, uJ, x\ = oo) and $i u (A;, u, X\ = oo). 
Those for $2s and $2u can be obtained immediately from ( [L.21 ). 

An important issue that arises at the outset of any calculation of universal scaling func- 
tions is that of proper choice of an ultraviolet cutoff. In the preceeding section, we chiefly 
used a Pauli-Villars cutoff to obtain the value of g c . However, we will see that g c does not 
explicitly show up in 1/N corrections and we are therefore free to choose the most conve- 
nient regularization scheme. All of the computations in this section were performed with 
two different cutoffs: 

(a) Lattice-cutoff: The ri£ field was placed on a square lattice of spacing 1/A, but no re- 
striction was placed on the allowed values of the Matsubara frequencies u n . The scaling 
functions were obtained in the limit A —>■ oo. By construction, in this cutoff scheme rela- 
tivistic invariance is violated at short scales and is present only in the long-distance theory 
at T = 0. Consequently, the T = spin-wave velocity will be renormalized from its bare 
value, and care will have to be taken to express the scaling functions in terms of the fully 
renormalized spin-wave velocity. 

(b) Relativistic, hard cutoff: The momenta k, and frequencies u n carried by the ng field were 
restricted to satisfy 



Unlike the previous cutoff, this scheme has full relativistic invariance, and there will be no 
renormalization of the bare T = spin-wave velocity. 

All of the remaining discussion in this section will use the second, relativistic cutoff scheme; 
as a result we will not have to consider explicitly the renormalization of the spin-wave 
velocity. We emphasize however that all of our numerical computations have been carried 
out with both schemes. While many of the intermediate results of the two schemes were 
different, the final results for the universal scaling functions were found to be identical. This 
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k 2 + oo 2 < A 2 . 



(4.1) 
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agreement provides strong support for the complete universality of the scaling functions, 
and makes it virtually certain that there are no numerical errors in our computations in the 
quantum critical region. 

The basic strategy for obtaining the scaling functions is straightforward. We evaluate 
Xs and Xu to order 1/N as outlined in Section [111 B , and then invert Eqns ( |1 . 1 aj , |TT 1 b| ) to 



express $i Sj « in terms of Xs,u- We will also need in this procedure the T = value of the 
ratio Nq / p s . The structure of the T = theory for g < g c is discussed in Appendix [D] where 
we found 



N$ gS* 



Ps N 



1 - V log 



(4.2) 



where the number rj = 8/(37r 2 iV) will become the critical exponent 7] defined in Section |. 

We will begin by noting some of the significant issues that arose in the evaluation of the 
results of Sec |H[B] in the quantum-critical region. We will then proceed to a statement of 



the results for the various scaling functions. 

The first step was to develop a fast computer program for the rapid evaluation of ie n ) , 
in ( [3.56D for mo = QT. Simple power-counting shows that II is convergent in the limit of 



the cutoff A — > oo. However, it is not clear a-priori that it is valid to take the A — ► oo limit 
at this early stage. The point is that subsequent integrals will involve values of H(q, ie n ) 
with q, e n itself of order A. However, we have shown by a detailed consideration of the 
relevant integrals, that these potentially dangerous contributions from n cancel out in the 
final results for all universal quantities. Thus we will fearlessly evaluate n in the limit of an 
infinite cutoff. 

The evaluation of n began with analytic determination of the integral over q\ in ( ^.56[ ). 
The summation over uo n was then performed by a direct numerical evaluation. Terms upto 
some large value of uo n were explicitly evaluated, and the remainder were summed by fit- 
ting to an inverse power series in upto order 1/uj®. A very similar procedure was 
used to evaluate dli/dirtQ. Finally the results were checked against the following computed 
asymptotic expressions: 
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n(g,ze n 



8(g 2 + 4)V2 



| (24-g 2 )9 3 T 3 l-6 7 | Q ( T" 



1 dU{q,ze n ) y/E q 2 + e 2 n 




8nQT (q 2 + el) 2 + AQ 2 T 2 e 2 n 



(4.3) 



where 



6 3 Je 



X 2 



7 = — / dx— = 2.32414354317. (4.4) 



e er 



Next, we evaluated the self-energy, S, and the constant K 3 defined in Q3.70Q . It is not 
difficult to show from (|4.3j) that both these quantities are logarithmically divergent in the 
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limit A — > oo. Further the coefficient of log(A) can be easily determined analytically. We 
numerically evaluated the integral over the momenta and the summation over the frequency 
for a fixed A and subtracted the known log(A) term. The remainder was found to be inde- 
pendent of A for large A, and was the required finite part of the result. These computations 
yielded the following catalog of useful results 



dY,(k = 0,iZJ= 0) 
dk 2 
T,(k = i®, iuj) 

dE(k = iO,iuJ = 0) 



i] log 



0.25266 



N 



—Q 2 rj\og 



A 



+ 



0.22616 



T J N 
0.69400 



dk 

E(fc,u7» 1) = (k 2 + uJ 2 ) 



77 log 



N 
A 



T 



-^\og(k 2 + u 2 ) + 



9lT 2 N 



K 3 (T,m = QT) 
T 



■— log(^) +0.17800. 

2tt 6 It; 



(4.5) 



where iuj) = Y,{Tk,iToJ) /T 2 . 

A little care is required in inferring the 1/N correction to the mass m at g — g c . The point 
is that the critical coupling g c itself has 1/N corrections, and in addition to the correction 
Sm 2 in (|3.60|) , there is an additional T-independent shift to itlq itself. Accounting for the 



1/N correction to the value of g c , we find the following result for m at g = g c 



m 



m 



K 2 (T, m ) 



(4.6) 



where K x , K 2 are to be evaluated using ( |3.60| ), ( |3.61| ), and ( |3.62| ) at m = 6T. The above 
result for m 2 differs from that in Section |IIIB| by the T-independent correction i^i(0,0) 
which in fact ensures that Sm 2 (T = 0) = 0, as should be the case at the gapless critical 
point. 

Now we need Ki(T, m = 0T). Simple power counting in (|3.62| ) shows that K\ is linearly 
divergent for large A. Moreover, the linear A term is ^(0,0), which from ([4.6|) , must be 
subtracted out. However, upon using the explicit results in ( |4.3| ) for II and dU/dm in ( ^.62| ), 



one finds that K\ is in fact only logarithmically divergent! The absence of any terms of order 
T/(q,e n ) 2 (which are allowed by naive power counting) in the asymptotic expansion of IT 
was crucial in obtaining this surprising result. The dangerous T/(q,e n ) 2 terms are, in fact, 
present at all values of m other than BT. Even for this special value of m , there is linear A 
contribution to K\ from q, e n of order A. However it was precisely these contributions that 
were dropped when n was evaluated in the limit of infinite cutoff. Since we are interested 
only in Ki(T,mo) — 7^1(0,0), we blithely neglect such contributions, and simply evaluate 
Ki as defined in ( 3.62|) using the values of n and dU/drriQ obtained above. The integral is 
only logarithmically divergent and can be evaluated in a manner similar to S and K 3 . At 
the end, we obtained from this the needed result for m: 
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™ 2 ^2 n2, M\ 0.21346 , AfTS 
— = 6 2 + V e 2 log ^- j + (4.7) 

We will now describe the universal scaling results obtained by combining the above with 
the results of Sec [IirB[ The factors of log (A) were found to cancel out of all universal 
quantities. 



A. Correlation length and structure factor 



For the correlation length we obtained 

1 



^XAoo) = 6 1 + 



0.2373 



(4.8) 



The residue of the structure factor in the vicinity of the pole in the complex k plane at £/£ is 
contained in the following result for the scaling function Hi of the structure factor (defined 

HD) 



in 



Ei(k, oo) 



0.4415N 



N )k 2 + Xl(oo) 
for k near iXAoo). 



(4.9) 



Our numerical results for the full scaling function Hi(fc, oo) for real k are presented in Fig. 
Analytic forms can be obtained in the limit of large and small k 



i (fc,oo) = 0.860818 + 



0.3697 
N 



F'f 0.864674-^ 



k < 1, 



(4.10) 



E x (k, oo) 



T(l- v /2) 2 ^k 



1-7? ' 



k > 1. 



(4.11) 



This last result for the behavior of Hi(&, oo) for large k required knowledge of the asymptotic 
properties of the scaling function $ is for the dynamic staggered susceptibility discussed in 
Section |1V (J| below; the constant Aq is given in ( |4.17| ). 

For experimental comparisons, it is convenient to express S(k) directly in terms of k£, 
where £ is the actual correlation length. From ( |L2jj ), ( 1.29|) , ( |4.10D and (|4.11| ), we obtain 



S(k) 



N*{Nk B T\ v V5 £ 



2vrp s 



2 (i + k 2 ey /2 



0.1925N 
N~) 



x < 



, 0.100N e i2 f . i 

Af / a/5 



(fee)" 

V5 



i + 



0.267 

"a - 



(4.12) 
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B. Uniform susceptibility 



For the scaling function O l5 (see ( |1.24j )), of the static uniform susceptibility we obtained 



oo 



^1 

— log 

7T 



V5 + l s 



0.6189 

N 



+ 



5xi 



(4.13) 



We have also performed Monte Carlo simulations of a lattice version of the 0(3) sigma model 
which are described in Appendix |Cj. This yielded Qi(oo) = 0.25 ± 0.04, in good agreement 
with the above result. 



C. Staggered Susceptibility 

We first describe some asymptotic limits of the two-parameter scaling function 
§i s (k,u, oo). For small k we have 

Re$^(M,oo) = 6 2 fl4 01830 



N 

T 2 ( 0.0001 \ -r- 
+ k (1 - -j^j ; k < 1. (4.14) 

Our numerical accuracy is not sufficient to rule out zero 1/N correction to the coefficient of 

— 2 

k ; the small value obtained for this correction appears to be accidental. At large arguments 
we found 



$J-(fc,i57,oo)= (k 2 +ZJ 2 ) 



l-^\ og (k 2 + uJ 2 ) 



+ '/ ( ( ^) + 3 



fe,w»l. (4.15) 



We expect that this logarithmic series can be exponentiated. Performing the exponentiation, 
followed by an analytic continuation to real frequencies we get finally 

$ ls (fc,a7,oo)= Aq ; £,^»1, (4.16) 

where the universal number Aq is given by 

A Q = 1-77 (log ( 4 - 17 ) 

to order 1/N. 

We turn finally to the determination of the scaling function $i s (&, uJ, oo) for arbitrary 
fc, aJ. As analytic evaluation is clearly impossible, we will have to resort to numerical 
computations. Moreover, as it is not easy to analytically continue numerical data, we will 
perform the numerical computations directly at real frequencies. There are some interesting 
issues which arise from the interplay between the analytic continuation, finite-temperature 
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effects, and ultraviolet divergences. These issues do not appear to have been discussed 
elsewhere before, and it appears worthwhile to present some details. 

Our strategy will be as follows: we will start with the key observation that to order 1/N 
ImS is free of ultraviolet divergences and that 



Im$ ls (k, ui, oo) = ImS(/c, a;) 



(4.18) 



Therefore Im^^ 1 can be obtained by a direct evaluation of Q4.18 ). All ultraviolet divergences 
and A dependences are in fact in the real part of self-energy. Next we note that ^~[} 
is analytic as a function of 57 in the upper-half frequency plane. However, we will find 
that because Im^^ 1 ~ 57 2 for large 57, its Kramers-Kronig transform is not convergent and 
cannot be directly used to obtain the real part. Instead, we will use the analytic information 
contained in the large 57 behavior in ( |4.15| ) to perform an appropriate subtraction, and take 
the Kramers-Kronig transform of the remainder. 

First some details on the evaluation of ImS. From the results for £ in Section |111 B , we 
obtain after and analytically continuation to real frequencies 



ImS(/c, uj) 



1 



4?r 2 A 



d 2 q 



dQlm 



o 



n(?,n). 



X- 



c q+k 
+ (1- 



-n|) 



(4.19) 



Here 



'k 2 + 9 2 ) 1//2 , oo > (ImS is an odd function of uj 
function, and there is no cutoff in the q integration. The propagator l/U(q,Q) is simply 



n(ej-) is the Bose 



1/n in rescaled variables. Thus from (|3.56| ) we get Imn 



Imn(g, H) 



^ 1 



d 2 p - 



e p+k e q 



167T 

n p+k~ n pl 5 (^ 



l C p+k 



(4.20) 



where Q > (Imn is an odd function of Q), and there is no cutoff in the p integration. 
The real part Ren can be obtained by a Kramers-Kronig transform of Imn (the frequency 
integral is free of ultraviolet divergences), which can then be used to determine Im(l/Tl). 
Note all of the integrals above defining ImE are pure numbers and amenable to a direct 
numerical evaluation, which we undertook. The presence of delta- functions in the integrand 
considerably speeded up the numerical computations. 

oo 



Next, we turn to Re$ ls . We deduce from ( [4.15|) that for momenta k fixed, but 



we must have 



Im$ ls 



- (u 2 logu 2 - A q 1 uj 2 - n x {k) \ogu 

+ /i 2 (F) + ..., 

Kl] _ /_ 2 /is 

— -sgnu; \u - Hx(k) + . . . 



(4.21) 
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where the functions ^i^(k) are unknown. We fit the numerically computed Im$ ls 1 to the 
above asymptotic form, and thence obtained H\(k). Then we defined the function P(k,UJ) 
by 

P(k,UJ) = Im^fcuJ) + ^sgn(aJ) (u 2 - ^(k)) . (4.22) 

The terms on the right-hand side have been chosen so that P falls off sufficiently fast at 
large UJ for its Kramers-Kronig transform to be well defined. Then, using the analyticity of 
$i s in the upper-half plane, we can conclude 

J-oo 71 11 — U 2 v 

-Aq'u 2 - ^(k) logul 2 ) + ^(k); (4.23) 



this determines Rec&^fc, to) upto the additive frequency-independent function fi2(k). Fi- 
nally, ^{k) was fixed by evaluating Re^^fc, To = 0) by an independent method: we deter- 
mined it directly from the expression ( |3.55| ) - unlike the computations just discussed, the 
frequency sums were evaluated by a direct summation along the imaginary frequency axis 
and a straightforward numerical evalutations of the relevant Feynman integrals. 

This completes our discussion of the derivation $i s (/c, u;, oo). The numerical computa- 
tions are summarized in Fig. [3]. 

One important feature of Fig. || which hasn't been discussed so far, is the behavior of 
Im$ ls at small frequencies. This can be determined directly from the expressions ( [4.1 9| ) and 
Q) for Im£. We found 

if 30 2 \ 

Im$ ls (fc = 0,57,oo) ~ — exp f "^J , ( 4 - 24 ) 

while 

-If 36 2 \ 
Im$ ls (|cJ| < k, oo) ~ ^exp f J . (4.25) 

These peculiar singularities are probably artifacts of large N expansion, because they arise 
from the strong contraints imposed by the delta functions in ( [4.19| , f4~20l) and the difficulty 



in satisfying them at small momenta; in other words energy-momentum conservation dras- 
tically reduces the phase space for emission/absorption of spin-waves with the spectrum 
ZJ = (k + 2 ) 1 / 2 . Actually, for self-consistent calculations, we have to include the damping 
of intermediate excitations, which lifts the restrictions imposed by the delta functions. This 
should probably give Im<3> ~ uJ for small UJ, as in naive expectations. Note that in an exactly 
solvable 1 + 1 dimensional model of a quantum phase transition, where analogous scaling 
functions can be computed, no such singularities appear0. 



D. Local Susceptibility and NMR Relaxation 

Having obtained scaling results for the staggered susceptibility, we can now easily obtain 
properties of the local susceptibility xl- The scaling function, Fi, for Ihixl was defined in 
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( |1.31| ). We determined Fi by performing the integration in ( |1.32|) numerically. Our result for 
N = 3 was shown in Fig. [5| From our numerical computation we find for small frequencies 



sgn{u) 



0.06, 



1-V . 



\uJ\ < 1. 



(4.26) 



The power of u at small uj is fixed by the requirement that lm(xL{^)) ~ ^ for small uj. For 
large uJ, we use the large k, uj result for $ is in ( |4.15| ), and ( |1.32| ) to deduce that 



F(uJ) 



\uj\ 1 - , i 
A, 



d 2 k 
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Aix 



sm 



tti]\ sgn(ZU) r\u\ 
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A Q sm(irr]/2) 
sgn(u) — v 



7T7//2 



my Jo (zj 2 -k y-v/2 
\m\ > l. 



(4.27) 



Thus Fi(m) tends to a universal constant for large m. 

From ( |1.36 ) and (^4.26| ) we find that the scaling function Ri(xi) for the NMR relaxation 
1.35|) satisfies 



m 



Ri(x\ = 00) 



0.06 

~n~' 



(4.28) 



E. Specific Heat 



We consider evaluation of the expression ( |3.71|) for the free energy density T at g = g c . 
We will of course only be interested in !F(T) — J^(0) which is finite as A — > 00. 
We will need the value of H(k, iuj) at T = 0. From the result ( ^4.3| ) 



H(k, %uj)\t= 



1 



q,LJ < A. 



=0 8{q 2 + uj 2 ) 1 ! 2 ' 
Using this result, and the representation ( |3.36| ) for g c we get from ( |3.71|) 
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Repeated applications of Poisson summation formula^ simplifies this result to 
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(4.29) 



(4.30) 



(4.31) 
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All the integrals and summations in this last expression can be shown to be finite in the 
limit A — > oo after using the asymptotic expansion of II in (|4.3| ). The absence of any terms 
of order T/(q,e n ) 2 in ( |4.3| ) is again quite crucial; the structure of the T 3 term in ( }4.3| ) is 
also such that all potentially dangerous log(A) terms cancel. The frequency summation and 
momentum integration in ( [4.3 1| ) were performed numerically and led to a result proportional 
to T 3 . We then evaluated the specific heat and obtained 

AN 

*!(oo) = — - 0.3344. (4.32) 





V. RENORMALIZED CLASSICAL REGION 

We now proceed to the calculation of the scaling functions in the region where the ground 
state is ordered (g < g c ) and the temperature satisfies NksT <C 2np s , i.e. X\ <C 1. Under 
these conditions, the Josephson correlation length £j ~ p' 1 is much smaller than hc/hsT. 
At the shortest scales, the antiferromagnet possesses D = 2 + 1 critical spin fluctuations 
which continue to be described by the scaling function $ ls in (|1.37|) . The crossover to the 
Goldstone region (see Fig. |2|) occurs at scales k£j ~ 1 (u£j/c ~ 1). In this regime, the 
dynamics is governed by rotationally-averaged spin-wave fluctuations about a Neel-ordered 
ground state. We will focus in this section on the second crossover (Fig. |^), which occurs at 
the correlation length £ (£ £j, Kc/(kBT)) to a fully disordered antiferromagnet. At scales 
larger than £, all correlations decay exponentially in space and the low-energy dynamics is 
purely relaxational. For a qualitative description, one can neglect quantum effects in the 
vicinity of this crossover, and study a simplified, purely classical version of the problem - 
the classical lattice rotor mo defi We will see, however, that for a complete quantitative 
description, quantum effects cannot be neglected at any k. 

A detailed study of the staggered spin correlations in the renormalized classical region 
was performed by Chakravarty et. a/.IIHl in the framework of the perturbative renormaliza- 
tion group approach for a classical rotor model. In this approach, one starts the description 
at relatively short spatial and time scales where there is a perfect short-range Neel order, 
and one can distinguish between the longitudinal and transverse susceptibilities. At these 
short scales, the thermal coupling constant g? = kBT/p s , which measures the strength of 
thermal fluctuations, is small (notice that in two dimensions, gr is a dimensionless quan- 
tity). One then performs RG calculations to see how g? grows at larger scales. The scale 
where the running coupling constant becomes of the order of unity is identified with the 
inverse correlation length £. At larger scales, perturbative approach is unapplicable but it 
is assume that £ is the only large scale in the problem and the behavior at k < £ is 
not very different from that at k = This assumption has been explicitly verified by the 
Bethe-ansatz solution of the 0(3) sigma model in 2 dimensions^!. 

The 1/N expansion, which we use here, attacks the same problem but from a different 
perspective. We start with the Green's function which satisfies the mean-field equation for 
the saddle point. The structure of the saddle-point equation in Section |II1 Aj shows that the 



symmetry remains unbroken at all finite T. The mean-field solution has a gap for quasipar- 
ticle excitations which we identified, at N = oo, with the inverse correlation length. In other 
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words, the correlation length is finite from the very beginning. An obvious consequence is 
that the spin susceptibilities are isotropic functions in the spin space: Xiji.1i w ) ^ This is 
indeed what we expect from the true scaling functions in 2D antiferromagnet at finite T (Eqn 
|1.10a|) . On the other hand, the temperature dependences of observables are not necessarily 
correctly reproduced by the infinite N theory. We will show below that the perturbative 
1/N expansion for x\ <C 1 is actually an expansion in 1/(N — 2) logxi. We assume that the 
logarithmic terms can be exponentiated; the 1/N expansion thus yields corrections in the 
form of extra powers of temperature: x\ 2 \ For iV = 3, the power 1/(N — 2) = 1, and 
there will be substantial changes in the temperature dependences of the observables, and in 
particular, of the correlation length. 

Most of our results for the staggered dynamic susceptibility agree with the results of 
the RG treatment for the classical lattice rotor model (some minor discrepancies with 
Chakravarty et. a/.0 are found however). We are also able to go beyond previous results and 
obtain explicit expressions for various pre-factors and scaling functions up to the two-loop 
level. For the uniform susceptibility, we have calculated the temperature dependence of \u 
in a quantum antiferromagnet and found a linear T dependence at low T, with a universal 
slope. We emphasize that the temperature dependence of Xu ls ^purely quantum effect. It 
was absent in previous studies of the classical lattice rotor modeltM which had (for N = 3) 
Xu a/3 = XJ-(<W — ( n a.np)) = (x± is the transverse susceptibility at T = 0). 

We now proceed to a more detailed discussion of 1/N expansion. As in the quantum- 
critical region, we will use the result Q4.2|) for the value of Nq/p s at T = 0, to provide a 
counterterm for the log A contributions to the universal function <3>i s . However, unlike the 
quantum-critical region, temperature no longer sets the overall scale for fluctuations, and 
we find it useful to introduce a function $i s (k,u) related to $i s by 

* n ^ ( hc V f Nk B T W it n 

<$> ls (k,oj)= — — — ®i s (k,oj). (5.1) 



k B T I \ 2np s 



In terms of this function 



N 2 ~ 

Xs (k,cu) = ^$ ls (k,Lu). (5.2) 

Ps 

Next, at X\ 1, the typical frequencies u ~ c^^ 1 are much smaller than ksT/h and 
equal-time structure factor S(q) is simply related to $i s (k,uj = 0): 

k TN 2 

= k B T Xs {k } u = 0) = — °-$is(k, u = 0). (5.3) 

Ps 



The result for $ ls to order 1/N follows from ( |3.10D and (| 
&u(k,u) = (l-rylog AA 



16p s 

(k 2 + u 2 + m 2 + ^{k,u)) , (5.4) 
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where self-energy and mass renormalization terms (m 2 = + 5m 2 ) are to be calculated as 



in Sec IIIB 



We will now describe, in outline, the computations to order 1/N for £, m, and the 
structure factor. We will then proceed, in the subsequent subsections, to present more 
precise and detailed results for these and other static and dynamic observables. 



At N = oo, we have from (|3.27|) 



m °"l^ exp riv^J' (5 ' 5) 

where p^ = °° is given by ( p.250 . We emphasize that m is expressed here in terms of the 
T = spin-stiffness constant p s computed at N = oo. There are however 1/N corrections 
at T = as well, and the fully renormalized p s indeed differs from Q3.25| ). When reexpressed 
in terms of the total p s at T = 0, mo by itself acquires a correction of the order of 1/N; this 
correction will be important later as a counterterm for the leading ultraviolet divergence of 
m = 

At N = oo, the value of the correlation length, £, is given simply by £ = l/m . Com- 
paring Eqn ( |3.27| ) with the results of previous perturbative approaches, we see that neither 
the numerical factor in the exponent nor the temperature dependence of the prefactor in £ 
agree with the results of the two-loop RG analysis of Chakravarty et. alM. As we already 
discussed above, this is not surprising because their analysis was done for the particular 
case of N = 3. The two- loop /3-function for arbitrary N was first calculated by Brezin and 
Zinn-JustinEl in a perturbative expansion about the ordered state, and their result for the 
correlation length is 




^ ^ (5 _ 6) 



At N = oo, this expression coincides with mo, as it should. Further, if we formally expand 
the r.h.s. of ( |5.6|) in 1/N, we find terms of the form (1/(N— 2)) log(fcsT/ (hcm )) and (1/(JV— 
2)) log (\og(k bT /homo)). We therefore anticipate that the same terms should appear in the 
1/N expansion. This is indeed what we found in our calculations, as we now demonstrate. 

We first observe that all \og(ksT /hcmo) terms in the 1/N expansion come from inte- 
gration over spatial scales which are much larger than the scale set by the temperature. 
Accordingly, performing the calculations with the logarithmic accuracy, we can restrict the 
evaluation of £ and Sm 2 to a classical lattice rotor model; this implies that we restrict the 
summation over uj to the contribution at u — only and set ksT/hc as the upper cutoff 
in the momentum integration. This substantially simplifies the calculations and we easily 
obtain from ( [3.56 ) 



u , m k B T q+^q 2 + Am' 

°) = / log y — . (5.7) 

irqjq 2 + 4mg Am o 



Substituting this expression into ( |3.55| ), performing the integration, and using (|5.5|) for m , 
we obtain 
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N 



log X k + . . . ; A* 



\og[k B T/(hcm )] 
log [V^ 2 + m 2 / m ] 



(5.8) 



where dots stand for the terms of higher order in 1/N, nonlogarithmic classical contributions, 
and for quantum contributions. Note that as written, ( |5.8| ) is valid only for k 3> mo; for 
k = O(m ), we have with the logarithmic accuracy A^ = A m = \og[k bT / (hcmo)]. We now 
substitute Q5.8|) into (|3.61|) and, using the smallness of hcm /kBT, obtain after some algebra 



5m 



ma 
N 



-4 log -^ + 3 log A r 

ncm 



+ 



(5.9) 



Next, we have to show that (i) the actual expansion holds in 1/(N — 2) rather than in 
1/N and (ii) the logarithmic series are geometrical and therefore can be exponentiated. In 
principle, to prove any of these statements, one has to examine the structure of the pertur- 
bative expansion up to infinite order in 1/N. This is a rather difficult problem to analyze 
and we will be content with demonstrating that both expectations are consistent with the 
perturbative results up to order 1/N 2 . Specifically, we computed logarithmic correc- 

tions to E fciQ . The computational procedure is tedious but straightforward. We followed the 
general procedure described by PolyakovS: we identified various regions of virtual momenta 
which contribute to o with logarithmic accuracy and, evaluating the integrals, obtained: 



k 2 + S fc n = k 1 



l + ^(l + | + ...)logA fc 



2N 2 



log 2 X k + 



k\\k) 1/{N - 2 \ 



(5.10) 



precisely as we anticipated. We didn't perform 1/7V 2 calculations for 5m 2 , but it is very 
likely that the expansion for the mass is similar to that for E^o- We assume that this is 
the case, and assemble the 1/N corrections to m 2 into logarithmic and double logarithmic 
series. This yields 

2 , (k B T\ 2 
log 




heme 



kjr 

he 



hcm 

Mr 



N-2 

2N/(N-2) 



(x m f {N - 2) 



(5.11) 



In writing the last line we also assumed that the value of mo in double logarithmic terms 
can be replaced by the total mass m. The verification of this assumption requires a compu- 
tation of finite contributions to order 1/A^ 2 which we didn't perform. We now assemble the 
contributions in fl5.4Q , ( |5.10|) and ( |5.11|) and obtain 



$ ls (h,uj = 0) oc 



yl/(N-2) 

k 2 + £- 2 



(5.12) 
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where with logarithmic accuracy 



he \ (k B T(N - 2) 



l/(7V-2) 



2np s 



exp 



27rp s 



(JV - 2)k B T / 



(5.13) 



This agrees with the two-loop RG result 
ture factor is 



. From ( |5.12| ) and (|5.3|), the equal-time struc- 



S(k) oc£ 



2 , k B T(N-2) \ ( N -^ N -V 
{ 2tt Ps ) 

(l + |log(l+(fcQ 2 )) 1 ^ 2 ) 

i + (kO 2 



(5.14) 



For N = 3, this coincides with the result of Chakravarty et. a/.@. 

An advantage of the 1/N expansion is that we can go further than ( p.l2| - |5TT4| ) and 
calculate not only logarithmic contributions in 1/N, but also the regular ones. For this 
type of calculation, the classical approximation is not sufficient and one has to consider the 
full quantum problem. The computations are lengthy but straightforward. We will skip 
the details here: interested readers can obtain a fuller description of the intermediate steps 
directly from the authors. We will list here only a catalog of the results similar to Eqn (|4.5|) 



m 



i-log^ + ±l gl g ML + l(31og2-l + C + 0.384i; 
N nemo N nemo N 



(5.15) 



m + E(fc = imo, 0) = Z m 



1 + 



dT,(k -> 0,0) 



1 - l loglog H; - 1 - 3841 



dk 2 

k 2 + E(fc>r 1 ,o) = zk 2 

dE(k = im ,0) _ x 
1 dk 2 " Z 



1 + A7 loglog ^ + I ' 3518 
1 - llog(l + ilog(l + (A;0 2 ; 



2 , , k B T 0.2385 

iv loglog ^ + — 



2 , , k B T 1.9561" 
N l ° glOg he^- + — 



1 + rj log 



NA 
16 p s 



1 , , k B T 0.9561 
— log log — 

N S & hcm N 



(5.16) 



Here C is the Euler constant C ~ 0.5772. As before, we use these results to evaluate 
universal functions for various observables. 
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A. Correlation length and equal-time structure factor 

We start with the scaling function for the correlation length. From (|5.15 , |5.16|) , we find 



he \ (k B T(N-2)\ 1/{N - 2) 



k B T J \ 2np 

(N - 2)k B T 



xexpf-— ], (5.17) 



where £o has a rather simple form 



&=1 _ (31o g 2- 1+ C) (51g) 

The same result was recently obtained in a different way by Hasenfratz and Niedermayer@. 
They deduced £ by comparing the free energy of the Heisenberg antiferromagnet in moderate 
magnetic fields with the Bethe-ansatz solution for the O(N) sigma model00. Moreover, 
they argued on the basis of the numerical analysis, that the 1/N result for £o is in fact the 
first term in the 1/N expansion for 

/ e x l/(7V-2) 

&=(§) r(l + l/(iV-2)), (5.19) 

where T(...) is the Gamma-function. The numerical evidence for ( |5.19| ) is rather convincing 
and we will use ( |5.19| ) for the experimental comparisons in Sec |V11| . 

Further, the equal-time structure factor is given by ( |5.3|) and using (|5.16|) we obtain 

S(k) = S(0) f(H), (5.20) 

where 

/ k B T\ f(N -2)k B T\ l/{N ~ 2) 



27T psj \ 2np s 

xffl^l, (5.21) 



S(0) = 2irN£ 



N 

and f(k£) is a universal scaling function introduced first by Chakravarty et. aiJll, for which 
we obtain 



N \ N J \ 2 

(5.22) 

At /c£ 3> 1, i.e., in the 'Goldstone' region, ( |5.20| ) reduces to 
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For N = 3 this agrees with the result obtained by Chakravarty et. al£3 (their definition for 
S(k) differs from ours by a factor of N). Clearly, Eqn ( 5.23|) is nothing but the rotationally 
averaged result for the Neel-ordered antiferromagnet. In the ordered state, there are N — 1 
gapless transverse spin waves which contribute to the low-energy part of S(k) and one 
longitudinal component of fluctuations which has a finite gap and does not contribute at 
low energies. Rotational averaging then gives a factor (N — 1)/N, as in ( |5.23|) . 

Though our results are very similar to the scaling theory^, we observe that the two limits 
in ( |5.22|) cannot be assembled into the single one-parameter interpolation form for f(k^) 
suggested by Tyc et. a/H. This is not surprising however because the one-parameter scaling 
function was introduced as a convenient, but approximate way to interpolate between fc£ ^> 1 
and where the behavior of S(k) is known from the hydrodynamic considerations. 

Further, we emphasize that p s and Nq in ( [5 . 2 1| ) are the fully renormalized spin- stiffness 
and sublattice magnetization at T = 0. Only in this case, all ultraviolet log A divergences in 
1/N corrections to £ and S(k) are canceled out. Finally, in obtaining the universal functions 
( |5.17| , |5.22|) , we actually didn't use the condition p s <C J. It thus appears that at least to first 
order in 1/N, the universal behavior holds for arbitrary p s , i.e., in the whole renormalized 
classical region. This is a remarkable property of the quantum sigma model, and the the 
arguments that the universality at all a may hold at arbitrary N were elegantly displayed 
in the analyses of Hasenfratz et. alEMIl In any event, this implies that our results for the 
renormalized classical region, which were obtained in a theory valid near g = g c , are in fact 
valid at small T for all g < g c . 



B. Uniform susceptibility 

We turn next to the calculation of the static uniform susceptibility. The expression for 
Xu valid at arbitrary x\ was given in ( |3.68| ) and ( |3.70|) . In the renormalized classical region, 
it is convenient to introduce $i„ = ksT $i u (0, 0, x\) such that 

X„(0,0) = (^) 2 <£ lu . (5.24) 
he 

We then use n mo ~ k B T '/homo, and obtain for $i u 

S lu = - K 3 - ^ (m 2 - ml), (5.25) 

where given in ( |3.13|) , is the contribution at iV = oo and the remaining terms are 1/N 
corrections. For $^ we have 

*S = -2W/0»^^(log^ + 1 ). (5.26, 



Using the expression for mo at small xi, we can rewrite the N = oo result for Xu as 

X s :-(^) 2 (lx± + ^k B T), (5.27) 
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where x± = PsC 2 is the transverse susceptibility. We see that in the limit of vanishing 
temperature, Xu * s precisely the rotationally averaged uniform susceptibility of the ordered 
antiferromagnet. This is likely to be the exact result, and we expect that all 1/N correc- 
tions at T = can be assembled into the renormalization of x±- On the other hand, the 
temperature dependence of Xu ls a purely quantum effect (it is entirely due to the second 
term in (|5.26|) ) and 1/N corrections to dx S u/dT are indeed possible. We also observe that 



the slope ofyj versus T at iV = oo is twice as large as in the mean-field Schwinger boson 
approacht3E3. This is not surprising however, because the mean-field Schwinger boson the- 
ory is the N = oo limit for a a- model of the A-component complex unit field defined on the 
(jpN-i ma nifold. This model is isomorphic to 0(3) sigma-model only at a particular value 
of N = 2 and there is no reason why the iV — > oo limits of the two theories should be the 
same. One of main technical points of this paper is that the 1/N corrections in the 0(N) 
theory are numerically quite small, making it a superior approach to make precise numerical 
predictions. 

The corrections to $i„ were computed in the same way as for the correlation length and 
equal-time structure factor. We skip the details of calculations and discuss only the results. 
As before, we found that all divergent contributions in the ultraviolet are canceled out when 
the the result is expressed in terms of the fully renormalized transverse susceptibility at 
T = 0. We have explicitly checked that there are no other corrections at T — > besides 
the renormalization of x±- Moreover, we did not find any logarithmic corrections to the 
temperature dependent part of Xu which might have change the power of the leading In- 
dependent correction. This result is likely to be valid to arbitrary order in 1/N although 
the proof is lacking. At the same time, we did find the finite correction to the $i n and our 
result for Xu valid to order 1/N is 



= m 2 (lx, + (^) k S) ■ (5.28) 



TIC 2 



The renormalization factor, (N — 2)/N, of the linear T term is likely to be correct to all 
orders in 1/N. The argument is that in the XY case (N = 2), the spin-wave analysis 
for Xu i s f ree f rom divergences^ and predicts a cubic, rather than linear dependence on 
T Xu = + 0(T 3 ). We also note that for N = 3, the renormalization factor is 1/3 and 

the slope of Xu * s therefore significantly reduced from the N = oo result and now differs 
substantially from the slope of Xu ^ n the quantum critical regime. We will use this later in 
Sec. |VII| for the interpretation of the experimental data. 



C. Dynamic staggered structure factor 

The calculations of the dynamic susceptibility proceed along the same lines as in Sec.jTV]. 
We use integral representation of the polarization operator continued to the real axis, and 
compute the retarded self-energy. The calculations are rather lengthy, so we skip the details 
and focus only on the results: details of the intermediate steps can be obtained directly from 
the authors. At k£ 3> 1, we obtained 

„ \ N-l (hc k ) 2 N 2 1 , 
Xs (k,u) = —— 1 k) ^ 2 , — ^ 5.29 
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where e k = (hc k ) 2 {k 2 + £" 



Ps/Xx- We have introduced here the k— dependent 



spin stiffness, p k s1 and spin susceptibility, x±> which are values of these observables at the 
momentum scale k. In the present 1/N expansion, they are given by 



Ps 



(N - 2)k B T 
2ir 



g + - log (1 + (fcf) 5 



and 



x± 



2 N — 

N + N 



,Ps) 



(5.30) 



(5.31) 



The dimensionless, numerical variable g was found, to first order in 1/N, to be g = 1 and 
independent of the ratio of u/hck as long as this ratio is less than 1. Eqn.( |5T3~0|) (with 
g = 1) coincides with the one- loop result for the running spin-stiffness inferred from the 
RG equation for the static coupling constant!!. It was suggested by Tyc et. alf^ that the 
two-loop corrections may lead to g ^ 1. In our approach, the same is likely to happen from 
higher order terms in the 1/N expansion; dependence of g on the ratio u/Tick is also possible. 
Furthermore, we have checked that Eqn ( |5.31| ) coincides with the perturbative result for the 
transverse susceptibility measured on the momentum scale k (and frequency scale ck). The 
first term in ( [5 .31 ) is the exact result at T = 0, which we already obtained in the previous 
subsection. The second term (which, we emphasize, is also a classical contribution) actually 
accounts for the difference between transverse and longitudinal susceptibility measured at 
finite momentum and frequency. In this situation, we probe the system at finite spatial and 
time scales where the system appears Neel ordered. The temperature dependence of this 
term for N = 3 is the same as in the scaling approach0 and, in fact, can be deduced directly 
from the diagrammatic expression of Xu m Sec. |111B| if k and u are both small but finite. 
The 1/N result for the damping rate 7^ ^ is given by 



Jk,< 



TT KCkk 

2 N — 2 



\N-2)k B r 



2irp\ 



k B r 

Tick 



(5.32) 



where (=1 for on-shell excitations) is a smooth function of the ratio oo/c k k This re- 
sult for the damping rate agrees with the lowest-order perturbative calculations by Tyc and 
Halperinlizl. They, however, have shown that the logarithmic dependence on the quasipar- 
ticle momentum k in ( |5.32| ) is actually an artifact of the Born approximation. Significant 
corrections to the self-energy term arise from the damping of intermediate excitations. Ne- 
glecting this damping is a legitimate approximation only if the damping rate is smaller than 
the energy of the incoming quasiparticle. Let us define the momentum scale for intermedi- 
ate states, q m , such that •y grn = e k . Clearly then, the lowest-order calculations are valid for 
q < q m , but damping of intermediate states must be included for q > q™. A simple estimate 
based on Eqn. ( f5.32j ) yields q m ~ k(p k s /T) 2 . A careful consideration!^ then demonstrated 
that q m has to be taken as the upper cutoff in the momentum integral leading to ( |5.32|) , and 
A;— dependent logarithm in ( |5.32[ ) has to be substituted by the self-consistent expression 
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The self-consistent result is then 



_ 7T hc k k _ f(N-2)k B T\ 2 
21 ° g (iV^T +r 



(5.34) 



where % >aJ (=1 for on-shell excitations) is a smooth function of the ratio u/ckk, which 
generally differs from 77, and T is a pure number of order unity. 

We return to consideration of the dynamic susceptibility Xs(k,(j) but now at momenta 
other than fc£ ^> 1. At kt; ~ 1 the naive 1/N expansion is not terribly useful because the 
scale q m ~ 0(£ _1 ), and the damping of intermediate excitations cannot be neglected at 
any momentum. We have therefore little to add here to the results of Reffl who used 
the dynamic scaling hypothesis^! to study this region. The content of this hypothesis is 
that there is only one large spatial scale in the problem, namely the correlation length, 
and the damping of excitations becomes comparable to the real part of the spectrum at 
k^ = 0(1). Indeed, if we use our result (|5.34 ) for 7^, which is strictly speaking valid only 



for fc£ S> 1, and extend it down to k£0(l), this is precisely what happens. The ratio 'jk^/tk, 
which behaves as (T / p h s ) 2 log p k s /T at large k£, becomes of the order of unity at the scale of 
k^ = 0(1). This happens because p k s in ( |5.30| ) is renormalized substantially downward at 
small k and eventually becomes ~ 

Notice also that x± tends to a finite value as k approaches zero. As a result, Ck decreases 



with k, and for kt; = 0(1), we have Ck ~ c ^ksT /2np s independent of N. This last result 
illustrates an important feature of the 1/N expansion in the renormalized-classical region. 
The A^-independent result for Ck is inconsistent with our earlier N = 00 analysis in Sec [111 A 



which is completely Lorentz-invariant and has a T-independent value of c. This apparent 
contradiction is a consequence of the non-commutativity of the T — *• and N — >• 00 limits. 
The key point is that the T = static uniform susceptibility in the 0(A) sigma model 
scales as 0(1/N) (Eqn ( |5.28| )) and therefore vanishes at A = 00. In view of this, one has 
to proceed to the next order in the expansion over 1/N (as we did) to check that there is 
indeed a breakdown of Lorentz-invariance at small but finite T. 



D. NMR relaxation 



The assumption about the functional form of Xs{k,uj) at small A; is a key ingredient 
which makes possible the calculation of the antiferromagnetic contribution to the spin- 



lattice relaxation rate. Using the definition of 1/T\ in ( 1.34 ) and performing the integration 
in (|5.29D , we obtain 



i-Xs 



A\Nl 
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1/2 
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(5.35) 



where A is a numerical factor whose calculation requires us to know the precise form of 
Xs(k,(j) at = 0(1). The functional form of 1/T\ S in ( |5.35|) is identical to that obtained 
by Chakravarty and OrbachEH on the basis of the scaling approach of Chakravarty et al. 
They also estimated the value of numerical factor to be AA5 2 « 0.61 (for iV = 3) by fitting 
the scaling forms of Chakravarty et a/t3 and Tyc et al0 to the numerical simulations on a 
classical lattice rotor model. 



E. Specific heat 



A simple inspection of the free energy ( |3.71|) in the renormalized classical regime shows 
that the dominant contribution to JF(T) — JF(0) comes from the region of magnon frequencies 



comparable with the temperature. In this region, mo in ( 3.71|) can be neglected compared 



to k 2 + u 2 , and we obtain 



IL(k, iu) = -j^r 



N k 2 + u 2 



(5.36) 



The free energy ([3.71 ) is then easily seen to be precisely the same as that of a gas of N — 1 
gapless Bose degrees of freedom. By definition, ( |1.10c| ), \?i(T — > 0) is a measure of the 
number of such modes in the ground state. We have therefore 



*i(a:i < 1) 



N-l. 



(5.37) 



It can be shown using the above results that x\ dependent corrections to are suppressed 
by a factor e~ 1//xi . 



VI. QUANTUM DISORDERED REGION 

The section will present computations of 1/N corrections to the staggered susceptibility. 
We will not compute 1/N corrections to the uniform susceptibility and the specific heat: 
both these quantities are suppressed by factors of e _A// ^ sT ^ at low temperatures, and are well 
described by the N = oo theory in Section. |111 A 2| -1/N corrections will lead to innocuous 



numerical prefactors. 

As in the previous sections, the computations will be carried out with a relativistic cutoff 
scheme which restricts lj 2 + q 2 < A 2 . We will begin with a description of results at T = 0, 
followed by a discussion of the finite T corrections. 



A. T = 

An immediate simplification here is that all correlators are completely relativistic. In 
fact, as the ultraviolet cutoff is also relativistic, this is also true at frequencies and momenta 
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of the order of A. All Green's functions are therefore dependent only upon a relativistic 
3-momentum Q, related to the usual momentum, q, and real frequency, u, by 



Q 2 = q 2 - {u + is) 2 . 



(6.1) 



Here e is a positive infinitesimal. 

An important feature of Xs(Q) has already been discussed in Section [TT"B : there is a 
perfect spin-1 quasiparticle pole at u = y/q 2 + A 2 . From (|3.54j) we deduce that the residue, 
A at this pole is (to order l/N) 



<9£(Q 2 



-A 2 



9Q 2 



(6.2) 



We turn therefore to an evaluation of S. This will be obtained from the results of 
Section |III B| evaluated at T = with mo = A. From ( |3.56|) it is easy to obtain the exact 
result for II 



n(Q) 



iVctan (£ 
4ttQ V2A 



(6.3) 



Next, (|3.55| ) can be reduced to a one- dimensional integral for £(Q) 

1 M P 2 dP 
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P 2 + A 2 . 

From this result, numerical and analytic manipulations show 



(6.4) 



<9£(Q 2 



-A 2 



9Q 2 



rj log 



0.4817408231574 



N 



£(Q 2 » A 2 ) = Q 2 



'A' 



(6.5) 



The first result can be combined with ( |4.2|) to obtain the result ( [1.22 ) for Zq, while the 
second is closely related to the fourth equation in ( |4.5|) . 



B. T > 



We present here results for the local susceptibility Imxii^) & t very small u. Clearly, 
because of the presence of gap, A, we have Imx L (u; < A)| T=0 = 0. Thus the entire contri- 
bution below comes from thermally excited quasiparticles. From the definition (|1.16|) of xl 



and ( |3.54j ) we have to order l/N 



IrnxL^uJ -> 0) 



gS 2 f d 2 k lmS(k,u) 



N J 4vr 2 (P + A 2 ) 2 



(6.6) 
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Using the integral representation of the polarization operator, it is not difficult to obtain 
that in the quantum disordered region (T <C A) 



JmE(k,u -> 0) 



2u f d 2 q e~ tk +" /T 



e k+g 



NT J 4tt 2 



Imn l (q,e k+q ), 



(6.7) 



where, as before e k = y k 2 + tuq. Note that e k and II now have to be evaluated at m = A. 

It now remains to evaluate Imll -1 . We will begin by considering Imll. From the gener- 
alization of ( |4.20| ) to the quantum disordered region we find the following important contri- 
butions to Imll for T < A 



Imll(g, e k+q ) = lmU(q, e k+q )\ T=0 + 



+ - 



57T 



2 n p n p+q r/ 

d p -o(e p ^ 
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e n — e k 



(6.8) 



where n|r =0 was obtained earlier in ( |6.3|) . In principle, the T dependent corrections to Rell 
can be obtained by a Hilbert transform of (|6.8j ); however, we found that these contributed 
subdominant corrections to Imx^/u; in the limit T< A. 

We can now compute by combining (|6.3|) , ( |6.6|) , ( |6.7|) and (|6.8|) . We omit the 

details and give the final result 



Imxiiu — > 0) 
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Jog^(A/r), 



(6.9) 



Note that the first term in the brackets is due to the temperature dependent part of Imll, 
while the logarithmical terms come from the real and imaginary parts of n|y =0 . This result 
can be combined with the definitions in Section [1 D 5| to obtain the small argument limit of 
the scaling function F 2 . 



VII. COMPARISON WITH EXPERIMENTS 

In this Section, we compare our theoretical results with the available experimental data 
for undoped, and weakly doped La 2 - x Sr x CuO i and the numerical results for 2D S = 1/2 
Heisenberg antiferromagnets on a square lattice. But first, let us briefly summarize our 
findings. 

We presented above the general forms for uniform and staggered susceptibilities in a two- 
dimensional quantum antiferromagnet which has p s J. The explicit crossover functions 
were calculated at iV = oo and the 1/N corrections were examined in the limiting cases of 
X\ ^> 1 and X\ 1, where X\ = NksT '/(2iip s ) is a parameter which governs the crossover 
between renormalized classical and quantum critical regions (for ii < 1, the system is in 
the renormalized classical region, while for x% ^> 1 it is in the quantum critical region). We 
found that for large values of x\, the perturbative expansion is regular in 1/N. Moreover, 
the corrections were numerically rather small; so we expect that for the physical case of 
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N = 3, the results obtained in the first order in 1/N are already quite close to the exact 
values of observables. On the other hand, at small x\, the 1/N expansion is logarithmically 
singular - it holds in logXi/iV and eventually changes the leading singularity in some of the 
scaling functions at x% — * 0; the final low-T behavior is the same as in the renormalized- 
classical scaling theory of Chakravarty et.al@. The crossover between small and large-xi 
behavior should occur at x\ around unity, though not necessarily at the same x\ for various 
observables. For p s <C J, x± ~ 1 is within the validity of the long-wavelength description, 
and we expect that there should be a temperature range where our formulas for the quantum 
critical region describe the experimental data better than the renormalized classical theory. 
Strictly speaking, the renormalized classical theory should be valid only for x\ <C 1 when 
log^i terms dominate regular perturbative corrections. 

We now proceed to describe the data. Let us first discuss undoped antiferromagnet. We 
know from elementary spin-wave analysis that zero-point fluctuations are not divergent in 
two dimensions and therefore the T = renormalization of spin-stiffness, spin-wave velocity 
and sublattice magnetization come primarily from the lattice scales, where we can rely on the 
results of spin-wave calculations. At present, two-loop spin-wave expressions are available^. 
For S = 1/2, they yield the values of p s , x± an d No. which are practically undistinguishable 
from the results obtained in numerical simulations^: 

p s = 0.181J; X ± = (^) 2 ^! N = 0.307, (7.1) 

where a is the lattice spacing. Hence 2np s « 1.13J and we therefore expect that quantum 
critical expressions should work for k^T > 0.35 J^OAJ. The expected crossover temperature 
is indeed not very small, but it is still significantly lower than J; in other words we may 
reasonably expect out long- wavelength description to continue to be valid for X\ > 1. Notice 
that our result for p s differs from p s « 0.15 used by Chakravarty et.a/.ll!. The reason is that 
Chakravarty et.al. extracted the value of p s from one- loop spin- wave results extended to 
S = 1/2, while our estimate of p s is based on a two-loop spin-wave calculations. 

We now discuss what happens at nonzero doping. First, a mean-field analysis based on 
Hubbard model predicts that antiferromagnet ism is destroyed at arbitrary small concentra- 
tion of holesll. However, more sophisticated considerations show@S that there is in fact no 
discontinuity in the immediate vicinity of half-filling, and one needs a finite, though small, 
concentration of holes to destroy antiferromagnetic ordering. There are several possible sce- 
narios of the doping-induced loss of Neel order: 

(i) There is a T = transition from the Neel state to a incommensurate magnetically or- 
dered stateiHi; or 

(m) The Neel is state is destroyed by quantum fluctuations and the system enters a quantum 
disordered spin-fluid with commensurate correlations. Only at a larger doping do incom- 
mensurate correlations appear. 

A recent self-consistent two-loop calculation!! on the Shraiman-Siggiail model displays both 
sequences of transitions depending on the strength of the coupling between fermions and 
n— field. 

The two scenarios differ primarily in their predictions for the behavior in the non Neel- 
ordered state; in both approaches, the spin-stiffness decreases with doping and vanishes at 
the critical point, while the spin- wave velocity remains finite at the transition. Thus, in any 



54 



event, the main effect of small doping is simply to decrease the bare spin-stiffness of the 
antiferromagnet. Moreover, the self-consistent analysis mentioned above0@ shows that the 
holes do not modify the critical properties of a Neel to quantum spin-fluid transition. Thus, 
the universal scaling functions computed in this paper can be used unchanged to describe 
this transition in doped antiferromagnets. 

A further, important, complication that has to be considered in realistic doped antifer- 
romagnets is the effect of randomness. All cuprate antiferromagnets have randomly placed 
dopant ions which will perturb the properties of the two-dimensional antiferromagnet. Ran- 
domness has been argued to be a relevant perturbation near the pure phase transition of 
the 0(3) sigma modelEl. It thus necessarily changes the universality class of the fixed point 
and also pushes the transition to smaller doping concentrations. In Appendix |A| we have 
presented a phenomenological discussion of the expected scaling properties of the random 
fixed point. 

However, we may conjecture (though we have no strong theoretical arguments for this) 
that the effects of randomness are important only at low temperatures in the immediate 
vicinity of the quantum transition. At higher T the dominant effect of doping is solely 
the change in p s and the properties of the antiferromagnet should be controlled by the 
pure fixed-point. We will soon see that neglecting randomness at moderate temperatures is 
consistent with the available experimental data in weakly doped La2- x Sr x CuO^. Of course, 
this simple approach will eventually fail at large doping. Crudely, we may take the largest 
possible x as one where p s would vanish in the absence of randomness. We will estimate 
from the data that this x to be somewhat larger than 0.04. 

An important consequence of the decrease of p s with doping is that x\ becomes larger 
at the same T. Hence the crossover between renormalized-classical and quantum-critical 
behavior occurs at lower T. Consequently, there should be a wider temperature range where 
our predictions for quantum-critical region should describe experiments better than the 
renormalized-classical theory. 

We now consider separately the experimental data for various observables. 



A. Uniform susceptibility 

We start with the uniform susceptibility, Xu- This quantity has no logarithmic correc- 
tions in the renormalized classical region, and is therefore an ideal candidate to test the 
predictions of the 1/N expansion. We consider first the numerical results for Xu ^ n the 
square-lattice 5 = 1/2 antiferromagnet. There have been high-temperature series expan- 
sions^, quantum Monte- Carlo! and finite-cluster calculations^ of Xu- Their results all 
show that Xu(T) obeys a Curie- Weiss law at high T, reaches a maximum at ksT ~ J 
and then falls down to a finite value at T = which is closeS to the rotationally aver- 
aged 1/S result (ha/gp B ) 2 x S u(T = 0) = (2/3) (0.514/8 J) « 0.04/ J (a is the lattice spac- 
ing). The data at low T are not accurate enough to make a reliable theoretical fit but at 
higher T (0.35 J < T < 0.55 J), both series expansions and Monte-Carlo calculations re- 
port a linear temperature dependence of Xu- The best fit to the Monte-Carlo data gives 
{fia/ gpB) 2 JXu = 0.037sci(l + 0.775/xi) (Fig. |7|). We compared this behavior with the 
theoretical prediction for the quantum critical region, which over the range of x\ values 
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used in the figure, is well approximated by the large-xi formula Eqn (|1.24j , |3 . 3 5| , [4 . 1 3|) : 
(ha/g[i B ) 2 Jx S u = 0.037xi(l + a/x{) where a = 0.8 + 0{l/N). (We do not know the l/N 
corrections to a; we only know that the l/N corrections modify the 1/xi term to 1/x 1 ^ .) 
The agreement between the slopes of the two results is remarkable. On the contrary, the 
slope for the renormalized classical region, Eqn (|5.28|) , is 0.014:ri in clear disagreement with 
the numerical data at T > 0.35J. 

We consider next measurements of X„ (T) in weakly doped La 2 -xSr x Cu0 4: . The interpre- 
tation of the experimental data requires cautioit^H because one has to subtract Van-Vleck, 
core and dia- and paramagnetic contributions from valence fermions from the measured 
Xu{T). Besides, as we said above, the effects of randomness are clearly important at low T. 
The subtraction of extra contributions neglecting the effects of randomness was first done by 
Johnson^ who actually estimated the strength of the Van-Vleck contribution by assuming 
that at zero doping, the susceptibility should be the same as in the Monte-Carlo studies of 
2D antiferromagnet. Nevertheless, his results for Xu clearly showllli! that at small doping 
concentrations, the susceptibility is linear in T with the slope of (ha/ 'gpB) 2 JXu vs - x i about 
0.043 which is close to our result (0.037). Unfortunately, we cannot definitely conclude from 
the experimental data whether the linear behavior the with universal slope stretches up to 
lower T as the doping increases, which has to be the case if p s decreases with the doping. 
We will observe this effect, however, in the data for the spin-lattice relaxation rate. Note 
also that at higher doping concentrations x\ ~ 0.1, the experimentally measured Xu(T) vs 
T flattens^. However, at such x±, the system is already in the metallic phase, where our 
approach clearly has to be modified. 



B. NMR relaxation rate 

The simple Mila-Rice-Shastry modeS for hyperfine coupling in La^CuO^ predicts that 
the hyperfine coupling constant for 63 Cu is finite at the antiferromagnetic ordering momen- 
tum (it, 7r), where the dynamic susceptibility is peaked. The region around (it, it) thus gives 
the dominant contribution to e3 Cu relaxation and we may use the long-wavelength theory for 
experimental comparisons. At the same time, the value of the hyperfine coupling constant in 
the sigma model approach cannot be inferred directly from the Knight shift measurements 
as in microscopic theories^!. Instead, we have to integrate out all intermediate scales in the 
microscopic model for the hyperfine interaction, to obtain the coupling constant between 
the nuclear spin and the unit vector field in the sigma model. This renormalization is not 
singular, however, and the fully renormalized coupling, which we label as A n , should not be 
very different from the microscopic one. 

At very low T (x\ small), the system is in the renormalized classical region. For iV = 3, 
our theoretical result Eqn. (|5.35|) is the same as in the hydrodynamic theory of Chakravarty 
et.alE 

1= /2 A i|A? 1 V /2 (72) 

T x V3 h he \2irpJ ' { ' ! 

Here A is the numerical factor which is difficult to calculate analytically. Chakravarty and 
Orbachli estimated it to be AiVg w 0.61 by fitting the scaling forms of Chakravarty et.ali^ 
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and Tyc et.alxS to the numerical simulations on a classical lattice rotor model. 

On the other hand, at higher temperatures {x\ > 1) we expect the quantum critical 
theory to work and 1/Ti should behave as in Eqn ( |1.35| , |1.36| ) 



Comparing these two forms for 1/Ti, we observe that the predicted spin-lattice relaxation 
rate rapidly decreases with the temperature at low T, passes through a minimum at x\ ~ 1 
and then slightly increases with T. In practice, 7] is very small for clean systems (77 ~ 0.028 
and therefore 1/Ti should be nearly independent on temperature in the quantum-critical 
region. 

We now turn to the data. The spin-lattice relaxation rate for nuclei coupled to the an- 
tiferromagnetic order parameter of 2D S = 1/2 antiferromagnet was numerically studied in 
high-T series expansions^ and finite cluster calculations^. In both cases, 1/Ti rapidly de- 
creases with increasing T at high temperatures (T > J/2) and becomes weakly temperature 
dependent around T ~ J/2. In finite-cluster calculations, the subsequent growth of 1/Ti 
at lower temperatures has also been observed. Clearly, this behavior is consistent with our 
theoretical observations. 

Experimental measurements of 1/Tjin undoped and weakly doped La2- x Sr x Cu04, have 
recently been performed by Imai et. alBQ in the temperature range 20 — 900K. For undoped 
system, the behavior above T N = 308K, but below 700K is well fitted by ( |7.2j ) with A n = 
(1.33 ± 10) 10~ 2 K and J = 1590 ± U0K (our estimate for A w is slightly different from 
theirs^ because we use the exact prefactor for the correlation length). However, at about 
650K, 1/Ti flattens and remains practically independent of temperature up to 900K - the 
largest temperature reported in by Imai et. a/.0. This is indeed what we expect from 1/Ti 
in the quantum-critical region. Furthermore, the experimental T range over which 1/Ti is 
nearly T independent increases with doping: at x = 0.04 it stretches nearly up to 500 K. 
We interpret this result as an evidence that p s indeed decreases with doping, thus pushing 
the system into larger x\ for the same T. 

For a quantitative comparison with the data, we need to know the value of i?i(oo). Direct 
1/N calculations give i?i(oo) = 0.06/iV, which is too small to account for the experimental 
result for 1/Ti. However, we already observed in Sec. |1V C| that Im$ ls calculated to leading 
order in 1/N has peculiar exponential singularities at small frequencies (Eqn. ( |4.24| , |4.25| )) 
which are probably artifacts of the large N expansion. These singularities substantially 
reduce the slope of Ti(uJ) at the smallest uJ (see Fig. ||). On the other hand, no such 
low-frequency suppression of Tifel was found in numerical studiesii and in the experiments 
on weakly doped La compounds^. In view of this, it appears reasonable to estimate the 
value of -Ri(oo) from our result for the scaling function F\ at slightly larger frequencies. 
Inspection of the numerical result for Ti (Fig. |5|) shows that Ti(cJ) is linear in uJ for ZJ 
between 0.5 and 1, and the slope yields i?i(oo) w 0.22. Substituting this result into ( \l.3\) 
and using the values of A n and J from the low-T (renormalized-classical) fit, we obtain 
1/Ti = (3.2 ± 0.5) x 10 3 sec _1 ; this is in a good agreement with the experimental result 




(7.3) 



1/Ti ss 2.7 x 10 3 



sec 



-1 
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C. Correlation length 



Detailed measurements of £(T) in pure La^CuO^ have been performedc3t3E3 at low T, 
where the system is in the renormalized classical region. Here Eqn. ([5.17 ) is valid and using 
(|7.1|) we obtain 

£(T)^0.50aexp{±^}. (7.4) 



Chakravarty et.alx^ used a different numerical pref actor in ( 5.17Q and a different value for 



the spin-stiffness constant. The combination of the two yielded nearly the same value of 
prefactor as in (|7.4j) , but the numerical factor in the exponent was slightly different (0.94 
instead of 1.13). This discrepancy is not crucial however, and both Eqn. ([7.4j) and the 
analogous expression of Chatoava/ty et alE fit the experimental data betweeTsSOA" and 
560K rather well. The value of J has been estimated by high-energy neutron scattering 
measurements of the spin- wave velocity" to be J ~ 1560K. Hence x\ ~ T/590K, and 
all of the experimentally accessible temperature range is within the renormalized classical 
region. Nevertheless, we estimated the value of £ at the highest experimentally accessible 
temperature 560K by using the N = oo expression for the crossover function X\ in ( |3.26| ); 



we obtained = 0.023V4™ 1 , which is not far from the experimental valueE3 of 0.03v4 _1 . 

At finite doping, we expect that the crossover between the two regimes will occur at 
smaller temperatures; quantum-critical behavior should therefore be observable at temper- 
atures above and around 500i^. Deep in the quantum-critical region, we expect that £ 
behaves as 

r 1 = 1-039 ( 1 - 1) , (7.5) 



hca \ X\ 

where 7 ~ 1. The 1/N corrections have been included in the slope but are not known 
for 7: they also change the subleading term to x\ . We fitted the data of Keimer et. 
a/.H at x = 0.04 by this formula and found satisfactory agreement with the data over 
the temperature range between 300i^ and 550K. The value of p s extracted from the fit: 
2np s ~ 150K 3007^, is still positive, but of course is much smaller than at zero doping. 
Note that Keimer et. a/S, used a phenomenological form for £ _1 (T) which combined the 
renormalized-classical result at zero- doping and temperature-independent correction due to 
finite doping; this form agreed well with the experimental data for doping concentrations 
x — -T- 0.04, but the theoretical arguments behind it are unclear to us. 

We also compared our results for the quantum-critical region with the numerical data 
for £ in a pure 2D S = 1/2 Heisenberg antiferromagnetilSff Numerical simulations were 
performed up to temperatures of about 4J; if quantum-critical behavior for £ is present in 
the Heisenberg antiferromagnet in some temperature range, it should have been detected. It 
turns out, however, that up to /c#T ~ J. numerical data are well fitted by the renormalized 
classical theory, although the best fitE3 gives the value for the prefactor which is nearly 
half of that in ( |7.4j ). On the other hand, for 0.4J < ksT < 0.6J (where 1/Ti levels off to a 
constant value), our pure quantum-critical result for £ is close to the numerical on MB. The 
interpretation of the numerical data at higher T requires cauton as these data at ksT > 0.6 J 
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are equally well fitted by the quantum-critical result oc x±(l — 7x1) where 7 is close to 
oneB However, theprefactor in the fit is nearly twice as that in (|7.5|) . We argue in a 
separate publication^ that these discrepancies in the fit to the quantum-critical theory are 
chiefly due to nonuniversal corrections ~ T/J which cannot be neglected above 0.6 J. 

Separate Monte-Carlo calculations of the correlation length precisely at p s = have been 
performed by Manousakis and SalvadoiS. They simulated the quantum 0(3) sigma model 
directly, rather than the spin-1/2 antiferromagnet. Their results yielded the value of the 
universal function Xi(oo) = 1.25 which is not far from the result of Chakravarty et. a/.0, 
or our result in ( J7.5[ ): Xl(oo) = 1.04. 



D. Equal-time structure factor 



There are, to our knowledge, only few data available on S(k). The behavior of S(k) vs k£ 
in a 2D antiferromagnet was studied by a Quantum Monte-Carlo by Makivic and JarrellE3. 
They fitted their data at k B T = 0.35J-^0.42J by the renormalized-classical scaling formula 
of Chakravarty ei.a/.0. We present another interpretation for the data. The key point is 
that at k B T = 0.42 J, the correlation length is about5.7a and the magnon energy is therefore 
hck/kBT ~ 0.7(A;£). We see that it becomes larger than k B T already at k£ ~ 1.5, and the 
use of the classical description at larger k is hardly justified. We rather have to use the 
full quantum expression for S(k) to fit the data. Clearly, at k B T ~ 0.4J, the Josephson 
correlation length is not very different from £, that is at Tick > k B T the system should be in 
the quantum-critical regime. However, we found above that the 1/N corrections to S(k) in 
this regime are very small, even for N = 3, and for experimental comparisons we may well use 
the scaling function for S(k) computed at N = 00. For k B T = 0.42J, we have the following 
theoretical prediction, valid in both the quantum-critical and renormalized-classical regions, 
from flTggD and ( pO)) 



_ cothfa Jl + (k£) 2 } 
S(k) = S -. V K J , (7.6) 

'l + (kO 2 



where a = hc/2k B T£ ~ 0.35 and S = (iV 2 / ' p s )(hc£/2a 2 ) w 2.48. We fitted the Monte-Carlo 
data by (|7.6|), using the overall factor S as the only adjustable parameter; we found very 
good agreement with the simulations not only for hck > k B T, but for all k (Fig. |||). However, 
the value of S in the fit is 5^ = 3.61, which is somewhat larger than our N = 00 result of 2.48. 
As yet, we have no explanation for the discrepancy, and therefore cannot judge from the data 
whether at intermediate k the system is in the quantum-critical or in the Goldstone regime 
(S(k) in the quantum-critical regime is very close to the N = 00 result, while in the Golstone 
regime, S(k) chiefly differs from its value at iV = 00 by the factor (N — 2)/N = 2/3). At the 
same time, the good agreement we found in the momentum dependence of S(k) is clearly 
consistent with our conjecture that at k B T ~ 0.4J, the antiferromagnet is very near the 
crossover between renormalized-classical and quantum-critical regimes. 

Experimental data for energy-integrated S(k) are available!! for pure La 2 Cu04. Clearly, 
the experimental temperature dependence for the correlation length was inferred from these 
data. The experiments!! were performed in the temperature range of T < 560K, where 
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x\ < 1. The momentum dependence of S(k) was reported to be well described by a simple 
Lorentzian S(k) = 5(0)/ (l + (£;£) 2 ), and the temperature dependence of 5(0) agreed with the 
renormalized-classical result 5(0) oc T 2 £ 2 . The absolute value of 5(0) was not determined 
in the experiments, so we cannot compare the experimental result at highest accessible 
temperature with the theoretical expression in the quantum-critical region, as we did for the 
correlation length. The data on 5(0) at finite doping have not been reported, to the best of 
our knowledge. 



E. Local susceptibility 



Finally, we consider the experimental data on the momentum-integrated dynamical 
susceptibility Imxii^) — I d 2 klmx s {k,uj). Extensive experimental measurements were 
doneilS for La^Sro^CuO^ in a frequency range u = 2 -4- 45mev and for temperatures 
10 < k^T < 500K. They showed that the experimental data at all values of to and fc^T, 
obeyed the following functional form to reasonable accuracy 



Imx L (u;) = I{\u\)F(u/T). (7.7) 



This form is in agreement with the theoretical scaling form ( |1.31| , |1.32|) for clean systems, 
or (| A 1 0| , |A19| ) from Appendix A for random systems; further, we arrive at the theoretical 
prediction that I ~ |c<j| m . The exponent /i = r\ > in clean systems, while we expect \i < in 
random systemsIR Experimentally, it was found that I was approximately u independent 
at the larger frequencies; this is consistent with our results for the clean antiferromagnet in 
which fi — rj is very small. At smaller frequencies, I showed significant frequency dependence 
which could be well fit by an exponent fi = —0.41 ± 0.05. Such a behavior is clear evidence 
of the importance of randomness at low frequencies and low temper atures0. Note also that 
the scaling plot to test the uj/T dependence of the function F shows rather good collapse of 
the data for a wide range of values of w/T@I§, thus giving strong evidence of the presence 
of quantum-criticality. The universal scaling functions at the pure and random fixed points 
probably have a rather similar shape, and thus the presence of disorder does not effect the 
scaling plot very much. It is only the exponent fi which is particularly sensitive to disorder. 

We emphasize, then, that this experimental data clearly shows the presence of quantum- 
criticality and suggests that the effects of randomness are important only at low energies; at 
high enough u or k^T one can suscessfully fit the data at x = 0.04 by the quantum-critical 
theory for a clean antiferromagnet. 



VIII. CONCLUSIONS 

We conclude the paper by recalling some highlights of our results. 

We have presented the general theory of clean, two-dimensional, quantum Heisenberg an- 
tiferromagnet s which are close to the zero-temperature quantum transition between ground 
states with and without long-range Neel order. While some of our discussion was more 
general, the bulk of our theory was restricted to antiferromagnets in which the Neel order 
is described by a 3-vector order parameter. For Neel-ordered states, 'nearly-critical' means 
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that the ground state spin-stiffness, p s , satisfies p s <C J, where J is the nearest-neighbor 
exchange constant, while 'nearly-critical' quantum-disordered ground states have a energy- 
gap, A, towards excitations with spin-1, which satisfies A< J. The allowed temperatures, 
T, are also smaller than J, but no restrictions are placed on the values of ksT j p s or ksT/A. 

Our results followed from some very general properties of the T = quantum fixed 
point separating the magnetically-ordered and quantum-disordered phases. These properties 
are expected to be valid in both undoped and doped antiferromagnets, though not in the 
presence of randomness. They are (i) the spin- wave velocity, c, is non-singular at the fixed- 
point (i.e., the dynamical critical exponent z = 1); for the case of a vector order parameter 
this implies that the critical field-theory has the Lorentz invariance of 2+1 dimensions; (ii) 
on the ordered side of the transition, there is a Josephson correlation length, £j, related to 
the T = spin stiffness p s by the hypothesis of 'two-scale factor' universality which implies 
that p s = hcT/^j, where the number T is dimensionless and universal; (Hi) turning on 
a finite temperature places the critical field theory in a slab geometry which is infinite in 
the two spatial directions, but of finite length, L T = hc/ftsT), in the imaginary time (r) 
direction. The consequences of a finite T can therefore be deduced by the principles of 
finite-size scaling. 

Under these circumstances, we showed that the wavevector /frequency-dependent uniform 
and staggered spin susceptibilities, and the specific heat, are completely universal functions 
of just three thermodynamic parameters. On the ordered side, these three parameters are p s , 
the T = spin- wave velocity c, and the ground state staggered moment Nq; previous works 
have noted the universal dependence of the susceptibilities on these three parameters only in 
the more restricted regime of ksT <C p s . On the disordered side the three thermodynamic 
parameters are A, c, and the spin-1 quasiparticle residue A. 

We have calculated the universal scaling functions by a 1/N expansion on the O(N) 
quantum non-linear sigma model, and by Monte Carlo simulations. For p s finite, these 
scaling functions demonstrate the crossover behavior between the renormalized classical 
regime, when thermal fluctuations are dominant, to the quantum-critical regime, where the 
dynamics is govern by the renormalization group flows near the T = quantum fixed point. 
For small ksT/p s , the T-dependence of our results was similar to those already obtained 
by Chakravarty et aS. For large ksT/p s , most of our results were new. We found that 
the crossover between the renormalized-classical and quantum-critical regimes occurs at 
X\ ~ 1, where x\ = NksT '/2irp s . In a square lattice, S = 1/2 Heisenberg antiferromagnet, 
2np s m 1.13J, and quantum-critical behavior therefore should be seen at ksT > 0.4J. 
We compared our quantum-critical results with a number of numerical simulations and 
experiments on undoped and lightly-doped La2-sSr$CuO± in the intermediate temperature 
range, and found good agreement with the data, particularly for uniform susceptibility, 
NMR relaxation rate and equal-time structure factor. It appears, therefore, that the use 
of a 'small' p s point-of-view is quite reasonable even for a pure square lattice, 5 = 1/2 
Heisenberg antiferromagnet - while ordered at T = 0, this system is evidently close to the 
point where long-range order vanishes. A small p s approach also appears to be appropriate 
for lightly-doped antiferromagnets. 



61 



ACKNOWLEDGMENTS 



This research was supported by NSF Grants No. DMR-8857228 and DMR-9224290, and 
by a fellowship from the A. P. Sloan Foundation. We are pleased to thank R. Birgeneau, 
S. Chakravarty, D. Fisher, E. Fradkin, P. Hasenfratz, T. Imai, B. Keimer, M. Makivic, E. 
Manousakis, A. Millis, H. Monien, A. Ramirez, N. Read, A. Sokol and R. Shankar for useful 
discussions and communications. 

APPENDIX A: RANDOM ANTIFERROMAGNETS AND ANISOTROPIC 

SCALING 

This appendix will present a brief discussion of the extension of the results of this paper 
to the case of random antiferromagnets. We will restrict our attention to systems in which 
the randomness preserves the Heisenberg spin symmetry i.e. antiferromagnets with random 
exchange constants Jy. 

Random impurities induce perturbations on the clean system which are uncorrelated 
in the spatial directions but fully correlated along the imaginary time direction. This has 
the immediate consequence of breaking the long-distance Lorentz invariance of the pure 
system. It therefore becomes necessary to allow for anisotropic scaling in space-time at the 
quantum fixed point. It is conventional to introduce a dynamic scaling exponent z such that 
characteristic frequencies, u, scale with the characteristic wavevector k as 

uj ~ k z (Al) 

in the critical region. 

It is important to distinguish to two distinct classes of magnetically ordered phases that 
can occur in random antiferromagnets: these are (i) Neel ordered and (ii) spin-glass ground 
states. In the first of these the ordering moment has a definite orientation on each of the 
lattice sites. In the spin-glass, there is a long-lived moment on each site, but its orientation is 
random. We will consider the properties of the transition of these two magnetically-ordered 
states to a spin-fluid in turn: 

1. Neel order 

We present here some phenomenological scaling ansatzes for the quantum phase tran- 
sition from a Neel ordered state to a spin-fluid in a random antiferromagnet. The scaling 
arguments are similar to those employed for the superfluid to bose-glass transition by Fisher 
et. a/@, although they did not discuss the issue of universal amplitudes. We will restrict 
our attention to the magnetically-ordered side g < g c . The magnetically disordered side is 
expected to be gaplessa and its properties will not be discussed here. We also note that 
little explicit reference to randomness will be made here, its main role being the introduc- 
tion of a z 7^ 1. The results should therefore also be applicable to other quantum phase 
transitions in clean systems which have z ^ 1 and are below their upper critical dimensions. 
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We also assume below that the antiferromagnet is below its upper critical dimension and 
hyperscaling is valid. 

As in clean systems, we expect that T = ordered state is characterized by a Josephson 
length scale, £j, separating the Goldstone and critical regions. As g approaches g c , this scale 
must diverge as 

ij~(9c-g)- u . (A2) 

At scales larger than £j, the system should exhibit conventional Goldstone fluctuations with 
a well defined spin- wave velocity, c. However, because z ^ 1, the spin- wave velocity should 
exhibit non-trivial critical behavior as g approaches g c : 

c~(g e -g) v <-'- 1) . (A3) 
The ground state is also characterized by an average ordered moment N which vanishes as 

No ~ (g c - gf, (A4) 

and a spin-stiffness p s which satisfies^ 

Ps~{9c-9Y id+z - 2 \ (A5) 
where d is the spatial dimensionality. The exponent identity 

2(3 = (d + z-2 + r])u (A6) 

generalizes ( |1.11| ) and will be useful to us below. 

In the main part of the paper we showed that the three properties Nq, p s , and c, of the 
T = state, completely determine the entire finite temperature form of the susceptibilities 
and the specific heat in clean antiferro magnets. We argue below that, remarkably, this 
continues to be true even in random antiferromagnets which have z ^ 1. The value of c 
helps determined the appropriate scaling between space and time even though c itself has 
non-trivial critical behavior. 

The following discussion will specialize explicitly to the case of d = 2, although the 
generalization to arbitrary d is quite straightforward. We will also assume that z < 2, 
otherwise, the universality in the spectrum near the critical point will be broken by the 
higher order analytic terms in an expansion in k. As before, we will use units in which 
H = ks = 1. Of course, we are no longer free to set c = 1! 

The application of finite-size scaling to quantum systems requires that one determine two 
length scales characterizing the effects of (i) the deviations from criticality, and (ii) the tem- 
perature, and write down universal functions of their ratio. The length scale characterizing 
deviations from criticality is clearly £j. In clean systems, the length scale £t characterizing 
the effect of a finite T was given by £r = c/T. Scaling suggests that in systems with z ^ 1 
we should have £t ~ T~ l l z . The scaling functions will now depend on the ratio ^j/^t- Using 
( |A~2| ) we see that this ratio is measured by the value of (g c — g) uz /T. However from ( |A5| ) this 
is the same (in d = 2) as p s /T. Moreover, since this ratio also has engineering dimension 
0, there are no non-universal metric factors that can appear. Thus as in clean systems, the 
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scaling functions will have a universal dependence upon p s /T. For similar reasons, they can 
also depend only on uj/T. 

It remains to consider the wavevector dependence of the scaling functions. By scaling, 
this must appear in the combination q£j> ~ q/T l l z . Let us now try and determine the 
metric factor in front of this combination. There are two basic rules: (i) the metric factor 
should involve combinations of observables whose scaling dimension is 0; and (ii) the entire 
combination which appears in the argument of the scaling function should have engineering 
dimension 0. Once these rules are satisfied, one is guaranteed, by the principles of scaling, 
that no non-universal pre-factors remain. A little experimentation shows that the following 
satisfies these criteria: 

Zt = ^(-) 11/Z . (A7) 



We now have enough information to use the same steps as were used in Section |T| and 
obtain universal scaling functions of the observables. We will omit most of the intermediate 
steps here and go directly to the results. 

Consider first the susceptibility measuring fluctuations of the order parameter, Xs- We 
find 



N 2 f c\ 2 ( T ^ 2+ ( v ~ 2 ^ z 




(A* 




where $i s is a completely universal function and there are no non-universal metric factors. 
We have chosen the prefactors to satisfy the convention that scaling functions should remain 
finite as g approaches g c . One can verify from ( |A3| ), ( |A5| ) and ( |A6| ) that the pre-f actor of 
the scaling function, and the coefficient of q/T l l z are non-singular as g approaches g c . The 
scaling results for all the observables dependent upon Xs can now be obtained in a manner 
similar to that used for clean antiferromagnets. We will display explicit expressions for only 
two of them: the correlation length £ satisfies 

/Z / rp \ 

(A9) 

while the local susceptibility Ithxl is given by 

lM = ^LY'' F (»L). (A10) 

Ps \PsJ \T p s J 

Again X\ and Fi are completely universal functions chosen such that Xi(oo) and i*i(uJ, oo) 
are finite. 

The properties of the uniform susceptibility and the specific heat follow from an under- 
standing of the hyperscaling properties of the free energy density, T. A simple generalization 
of the arguments of Privman and Fisher^ to anisotropic systems yields (in d — 2) 
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J r = J r o + T^ 2 v(yj, (AH) 

where Tq is the ground state energy, and (p is a universal function. As in Section |I| the 
following results for static uniform susceptibility Xu an d the specific heat Cy now follow 

Xu = (9Vb 



'V 



T 









LX 


-2+2/z 


PsJ 






(A12) 



where are universal functions, with Oi(oo) , \&i(oo) finite. Note, in particular, that 

the Wilson ratio, W, (Eqn. ( |1.25| )) continues to remain a fully universal function of T/p s , 
and has a universal value at g — g c . The universality of the Wilson ratio, even in the 
presence of anomalous powers in Xu an d Cy, was also noted recently (using very different 
arguments) in the boundary critical theory of overscreened Kondo fixed points 



100 



2. Spin-glass order 

We now consider the properties of a quantum transition from a spin-glass ground state to 
a spin fluid. Clearly many of the properties discussed in the previous section are special to 
Neel ordered states and do not generalize. However measurements which are spatially local 
do have similar critical properties. As we do not wish to discuss the nature of the spin-glass 
state itself, we will restrict ourselves here to behavior in the quantum-critical region where 
T^(g c - 9 r. 

A very useful measure of the local spin correlations is provided by the local spin suscep- 
tibility, Xl{w) defined in ( |1.16|) . Along the imaginary frequency axis, this local susceptibility 
is given at the Matsubara frequencies u n by 



n) = / dre^C[r) 
Jo 



C(r) = <Si(0).Si(r)>, (A13) 

where the bar represents an average over all the sites i. The function C(r) can be used to 
distinguish the spin-glass and spin-fluid states. In the spin-fluid state C(r) will decay to 



zero for large r, will in the T = spin- glass phaset--l we must have 



lim C{t) = q EA > 0, (A14) 

T — >00 

with qEA the Edwards-Anderson order-parameter. It is conventional to define the order- 
parameter exponent (3 by the behavior of qEA as g approaches g c : 

qea ~ (g c - gf- (A15) 

The value of qEA thus fixes the behavior of C(r) at infinite time. We can obtain the behavior 
of C(t) for finite times r by a simple application of the dynamic scaling hypothesis. We 
expect at T = that 
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C(t) = (g c - gfhMo ~ 9c\n ! T = 0, (A16) 

where hi is a scaling function and z is the dynamic scaling exponent. Clearly we must have 
hi(x — > oo) be a finite non-singular constant to satisfy ( |A14[ ) and ( |A15| ). For r <C (# c 



standard critical phenomena lore requires that C(r) become independent of g c — g- This only 
possible if h±(x) ~ x~^l zv for small x. Putting this together with ( |A16j ) we can determine 
the behavior of C(r) at T = and g = g c 

C(r)~^y; T = 0, T<L{g c -g)- w . (A17) 



We can now use finite-size scaling to determine the behavior of C(r) at finite T 

C ^t) = ^jh 2 (Tr) ; T-\ r « (g c - g)-", (A18) 

where h% is yet another scaling function. Finally we use ( |A13| ), to determine xii^n)- In 
general we may find that the fourier transform is dominated by non-universal contributions 
at short times. However, upon analytically continuing to real frequencies, we expect all 
these non-universal contributions to affect only the real part of Xl, while the imaginary 
part is dominated by the universal long-time behavior. We therefore obtain the following 
generalization of ( |A10|) to systems with spin-glass order 



where 



Im. U M = kf/-j(£) : :/U< »(,,,.- ,,r. (A19) 



fi=-l + —, (A20) 
zv 



and F\ is a universal scaling function with a single non-universal overall scale. For small 
argument, we have F\(p) ~ sgn(u;)|u;| 1_/t , while F\{o — > oo) is finite and non-singular. 

In addition to the local susceptibility, the thermodynamic properties of the quantum 
phase transition which can be deduced from the hyperscaling properties of the free energy, 
are very similar in the spin-glass and the Neel ordered systems. In particular, we expect 
the following temperature dependence of \u an d Cv in the quantum-critical region (T ^> 
(g c — g) zu of the spin-glass to spin-fluid transition: 

V St rv T~ 1+d / Z 

C V ~ T d/z . (A21) 
It follows then that the Wilson ratio W (Eqn. ( |1.25| )) is a universal number at g — g c . 



APPENDIX B: BERRY PHASES AND DANGEROUSLY IRRELEVANT 

COUPLINGS 

We have assumed in this paper that the 0(3) sigma model is sufficient to determine 
the quantum-critical scaling functions of quantum antiferromagnets. A key step in this as- 
sumption is that the Berry phases present in the antiferromagnet can be neglected. This 
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assumption is based on the following circumstantial evidence. After some rather involved 
calculations, which have been discussed at length elsewhereElS, it was shown that the pri- 
mary effect of the Berry phases was to induce spin-Peierls ordering in the quantum-disordered 
phase of non-even-integer-spin antiferromagnets. This was shown in the context of large M 
calculations for SU(M) antiferromagnets. Further, a key feature of this calculation was the 
appearance of two well-separated length scales which characterized the fully-gapped spin- 
fluid phase. The first of these scales was the two-spin correlation length £ which determined 
the exponential decay of the equal-time spin-spin correlation function. The second was ^sp 
the length at which Berry phases first became effective in inducing spin-Peierls ordering. In 
the large M limit these two were found to be related by 



where Qi = 0.062296 + 0(1/M). It is clear that for sufficiently large M (which is the only 
region in which we know how to perform these calculations) we have £sp 3> £• 

It was then pointed out to us by Daniel Fisher^l that the appearance of two length-scales 
at a second-order phase transition, one of which is a power of the other, is a characteristic 
property of systems with dangerously irrelevant coupings. A dangerously irrelevant coupling 
is defined as one which is irrelevant at the critical fixed point separating the two phases, 
but is relevant at the fixed point which controls the nature of the phase one is studying. 
Crudely speaking, the coupling decays to a very small value at the first length scale while 
the system is controlled by the critical fixed point, but grows again to a value of order unity 
at the second, larger, length scale. 

For the antiferromagnet, our assumptions then, are the following. We assume that the 
dangerous-irrelevancy of the Berry phase effects, found in the SU(M) models, survives in the 
0(3) sigma model. The coupling to the Berry phase terms in the action decays to a negligibly 
small value after renormalizing out to a scale of order max(c/T, c/A). At scales larger than 
this, the coupling grows again, as the renormalization group flows are now in the vicinity 
of the strong-coupling fixed point controlling the quantum-disordered phase. However, this 
coupling significantly modifies only those correlation functions which are directly sensitive 
to the presence of spin-Peierls ordering. For all other spin-correlations (which includes all we 
have considered in this paper) the Berry phases can be neglected in determining the leading 
quantum-critical behavior. 

To clarify this issue, we now present a pedagogical discussion of a simple statistical me- 
chanical model with a dangerously irrelevant coupling. Unlike the quantum antiferromagnet, 
the overall structure of the renormalization group flows is well- understood in this model. 

We consider the finite temperature properties of a classical, XY model on a cubic lattice. 
At each site we introduce a four- fold anisotropy field /i 4 , which we will find is dangerously 
irrelevant. The model is described by the partition function Z 




Z 





(B2) 
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where the sites i,j lie on a 3-dimensional cubic lattice. This model will have a phase 
transition at some T = T c from a high temperature paramagnetic phase to a low temperature 
ordered phase. It is well known that the four-fold anisotropy h$ is irrelevant near T = T c 
and the phase transition is therefore in the universality class of the pure three dimensional 
XY model. However, it is also clear that the field surely cannot be neglected in the 
ordered phase. It breaks the 0(2) symmetry of the XY model, and must therefore destroy 
the Goldstone modes. Further, the common mean orientation at each site must be one of 
9 = 0, 7r/2, 7r, 37r/2 and cannot be arbitrary as in the XY model. The apparently conflicting 
properties of the critical point and the ordered phase are reconciled by the concept of a 
dangerously irrelevant coupling. 

Let us examine the structure of the renormalization group flows of this model for T close 
to T c : we measure the deviation from criticality by the reduced temperature variable 

T - T 

t =h jr~- ( B 3) 

A schematic of the renormalization group flows projected onto the T, H& plane are shown in 
Fig. p. The critical fixed point is at t = and = 0. The initial growth of t away from 
this fixed point is given by 

t{£) = te e/u } (B4) 

where < t <C 1, e e is the length rescaling factor, and v is the usual thermal critical expo- 
nent. The field h± is irrelevant at this critical point and will therefore decay exponentially 

h 4 {£) = /i 4 e"^, (B5) 

where h± is of order unity, to is the crossover exponent associated with h 4 at the critical 
fixed point. The system will emerge from the critical region at the Josephson length scale 
£j where t(£ = £{) « 1. From (|B3|) we see that 

6 = tfi = r v . (B6) 

For £ > £*, the flow of h\ will now be controlled by the T = fixed point. Let us assume 
that hi is relevant at the T = fixed point with eigenvalue <fi > (the value of will be 
determined later). Then we have 

h 4 (£ > £\) w h 4 (£\)e^-^ 

= h A e-M v+ M* + * e . (B7) 

Thus hi will return to a value of order unity when the argument of the exponent is zero. 
This defines a second length scale £ 4 = e e z where 

(B8) 

The effects of the h^ field thus become important at length scales of order £ 4 : this is the 
scale at which the Goldstone modes are destroyed, and the condensate gets locked at one of 
9 = 0, 7r/2, 7r, 37r/2. The scale £4 is related to £j by 
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£ 4 = £l+«V(<H ( B9) 

Additional insight can be gained by considering the scaling form for the transverse sus- 
ceptibility near the transition. Assume the condensate points at 9 = 0. Then the transverse 
susceptibility satisfies the scaling form 

IWI 2 ) = —eMk^cth^), (BIO) 

Ps 

where k is the wavevector, ip is a universal function, and C\ is the only non-universal metric 
factor. For most values of the term proportional to can be treated as a small per- 
turbation which makes a subdominant correction to the leading critical behavior. Only at 
extremely small values of k£j does the term become important. Matching to the expected 
form of the incipient Goldstone modes, we should have 

<p{k,Ti) = =2^= ] £<1- (Bll) 

k +h 

— — —1/2 

Thus the h term is significant for all k < h or for 

AT 1 > £ 4 ~ ^ ,{2v) . (B12) 



Comparing with flB~9f) we see that = 2. The key point, of course, is that £4 3> £j. In 
particular, the crossover from critical to Goldstone fluctuations occurs at a scale of order £j 
and is described by the scaling function of the pure XY model ip(k, 0). Only at much larger 
scales does it become necessary to include the effects of the field. 



APPENDIX C: MONTE CARLO EVALUATION OF QUANTUM-CRITICAL 

UNIFORM SUSCEPTIBILITY 

It was shown in Section [| that the high temperature behavior of the uniform, static spin 
susceptibility is given by (see ( |1.24|) 

xt(T)={^) 2 k B TQ 1 (oo). (CI) 

The universal number ^(00) has been evaluated so far in a 1/N expansion with the result 
( 4.1ij| ). In this appendix we will describe a determination of fii(oo) at iV = 3 by Monte 
Carlo simulations. 

The quantum 0(3) non-linear sigma model is expected to be in the same universality 
class as the classical, Heisenberg ferromagnet on a cubic lattice. Our simulations were 
therefore performed at the critical point of this latter model. We considered the ensemble 
defined by the following partition function 

z = jY[ds ie - n , 
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H = —K 22 Si ■ Sj, 

<ij> 



(C2) 



where i,j extend over the sites of a cubic lattice, and Sj = (S X) i, S y ^, S Z) i) is a 3-component 
vector of unit length. We used a lattice with L x L x L T sites, with periodic boundary 
conditions in all three directions. The Wolff single-cluster algorithmic was used to sample 
the states. This simulation was carried out at the critical value K = K c at which this model 
has a second-order phase transition. The value of K c is known very accurately from recent 
high precision Monte Carlo simulations 



i. 



K r = 0.6930. 



(C3) 



It has been argued that Xu * s related to the stiffness, p T , of this system to twists along 
the r direction. Upon examining the response of 7i to a field that generates rotations in the 
x — y plane, we obtain the following expression for p T 



1 



L 2 U 



K {S Xi iS Xj i+f + Sy^Sy^ 



l+T) 



l+T) 



(C4) 



where the expectation value is to be evaluated in the ensemble defined by Z. Finally, the 
universal number fii(oo) is defined by 



fii(oo) 



lim 



lim L tPt \ k=k 



(C5) 



It is crucial that the L — ► oo limit be taken first, to model a quantum system which is 
infinite in the spatial directions. The subsequent L T — ► oo places the quantum system at 
zero temperature. 

The results of our simulations are contained in Table |. Three independent simulations 
of 70000, 70000, and 210000, flips per spins were performed in systems upto L = 30 and 
L T = 10. A polynomial extrapolation to L = oo at L T fixed yielded the results shown in 
the last column of Table [[]. Finally, a second polynomial extrapolation to L T = oo was 
performed to yield the following estimate for fii(oo) 



fti(oo) = 0.25 ±0.04. 



(C6) 



APPENDIX D: COMPUTATIONS IN THE NEEL STATE AT T = 

In this Appendix, we derive the expressions for the T = sublattice magnetization A^, 
and the spin-stiffness p s to order 1/JV. These results will be used in the derivations of 
the universal scaling functions for the uniform and staggered susceptibilities in both the 
quantum-critical and renormalized classical regions. As in Sec |T|, our starting point is 
the functional integral for the O(N) sigma model. At zero temperature, the spin rotation 



70 



symmetry is broken and the perturbative 1/N expansion has to be modified to account 
for the nonzero expectation value of the order parameter. This 1/N expansion has been 
developed by Brezin and Z inn- Just inE! and we use some of their results. We first represent 
the unit vector field as 



n 



0o + 7T, 



(Dl) 



where (<7o) is finite and <?o-7r = 0. The sublattice magnetization N in spin-S antiferromagnet 
is expressed as iVo = S(ao), where S = SZs- The renormalization factor Z$ accounts for 
the order parameter fluctuations at short (lattice) scales which have to be integrated out 
in the mapping to sigma model from the original spin Hamiltonian on the lattice. Upon 
substituting (p 1|) into (|3.1|), the functional integral becomes 



Z = JvaoVm 5(a 2 + 7if-l) 



exp 



dr 



(V r a ) 2 + (V r vr,) 2 + ( (d T a ) 2 + (d T m 



(D2) 



where the index I now runs from 1 to N — 1, p s is the bare spin-stiffness and Co is the bare 
spin-wave velocity. All of the discussion in this Appendix will use the relativistic cutoff 
scheme when the momenta and frequency satisfy k 2 + uj 2 < A 2 . In this situation, the full 
relativistic invariance is preserved at each order in the perturbation theory and we will not 
have to consider explicitly the renormalization of the spin-wave velocity. To simplify the 
presentation, below we use the units where h = c = 1. As in Sec fTJ, we introduce a 
Lagrange multiplier A into the functional integral to impose the constraint, and integrate 
over 7T/. This gives 



VBqVX exp 



d r l dr 



(<Vo) 2 + K<% - !) 



N-l 




2 l0g 


1-31 + All) 



(D3) 



Here d 2 = (V r ) 2 + (d T ) 2 . The saddle-point point equation is easily obtained by taking a 
variational derivative over A and neglecting fluctuations in cr - This yields 



(N-l) 
N 



d 2 kdu 



G {k,iu) = 1 - (a )'' 



(D4) 



Here g = N/ p° s is the coupling constant and Go(k, iuj) = ik 2 + u 2 )^ 1 is the zero-temperature 
propagator of the 717 field at iV = 00. The 1/N corrections to ( [D4] ) are calculated as described 
earlieir3, with the modification that we have to consider the fluctuations of <To around its 
mean value. We find the T = analog to the polarization operator 



n*(£;, iuj) = Il(k, iuj) + - (cr ) 2 G Q (k, iuj) 

9 



(D5) 



and the correlator of the a field: 
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((a (k,ico)) ) — (a ) S(u})5 (k) 



G Q (k,iu) ( 2 2 G (k,iu)\ 

1 - - (^o) fF-r; • \ • ( D6 ) 



P° s \ 9 U.*(k,ico 
The condition a$ + 7if = 1, then yields 



l-(<xo>=£ / -— T G(fc,z^ 



(2vr) 3 

_ 2(g ) 2 r g? 2 £;gL; Gl(k,iuj) 

N J (2tt) 3 n*(ife,zw)" 1 J 

Here G is related to Go in the usual way: 

G-\k, iuj) = Go\k, iuj) + E(ife, iw), (D8) 



Below, we will express A" directly in terms of the fully renormalized spin-stiffness. It is 
instructive, however, to compute the critical exponent for N directly from ( P7|) . For this 
we observe that the first term in the r.h.s. in ( |D7| ) is simply a constant so that with the 
accuracy to l/N, 

U\* ( 9c-9 \ ( A | 2 f d 2 k du G 2 (k,iuj) \ 

{ao) = l^rj [ 1+ nJ wn^yj' (D10) 

where g c is the non-universal critical coupling. The value of the integral depends on the 
precise form of the polarization operator at k,u ~ A, which can be very complicated. 
However, for the logarithmic contribution in (|D10|) , we only need to know the form of 
U*(k,iu)) for momentum and frequency well below the upper cutoff. For such k and u the 
evaluation of H(k, iuj) at T = is straightforward and we obtain 

U(k,iu) = - 7 =±==. (Dll) 



iVk 



10 



We then use ( |D5|) for n* and evaluate the integral in (piO|) with the logarithmic accuracy. 
We obtain 

<->-(^)( 1 + ^^)- 

The correction can be exponentiated in the usual way and we get 

<a > 2 =(^j . (D13) 

where 
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2,3=1 (D14) 

Our next goal is to express N in terms of the fully renormalized spin stiffness at T = 0. 
At TV = oo, we have from ( p . 2 5|) 




(D15) 

Now we have to compute 1/N corrections to ( |13 1 5| ) . In principle, it is possible to evaluate p s 
directly by calculating the response to the twist of the order parameter in the momentum 
space. In practice, however, it is more convenient to calculate the static susceptibility x_l> 
which measures the response to the twist in r-direction; the value of p s then can be obtained 
using Lorentz invariance: p s = Cg • x±- 

The calculations of x± to order 1/N were described in Sec [Tl [T3[ As before, they have 



to be modified to account for nonzero order parameter. For N = oo, we use ( pi|) for the 
n— field, substitute it into the bubble diagram in Fig. [], and after simple manipulations 
obtain at T = 0, 

xl=°° = 4 (^o) 2 - (D16) 

4 



The calculation of 1/N terms proceeds along the same lines as in Sec [II [B[ We substitute 



(|DT|), into the 1/N diagrams in Fig. |6| and using ( |L)6| ) obtain after some algebra, 

2 

N 



X± = xl = °° (l-4^=o), (D17) 



where 



'T=0 



d 2 kdw G 2 (k,iuj) (k 2 -3m 2 ) 

~{2^f U*{k,tu) {k 2 + u 2 ) ' ( } 



We emphasize that the calculations at T = are much simpler than that at finite T because 
in fact we have to keep only the terms ~ (°"o); & H other contributions give zero after 
integration over intermediate frequency. This indeed is clearly seen from the expression for 
X± at finite T (Eqn ( |3.68j )) where each term contains derivatives of the Bose functions. 
Finally, from (pi3| ) and ( |D14|) , and the definition of Nq, we obtain: 

"S-^'l-iUwV (D19) 



Ps N \ N 



Clearly, the value of the integral in ( |D18| ) depends on the form of U*(k,iu) near the 
upper cutoff and the result for N^/ p s is therefore model dependent. However, we explicitly 
showed in Sec [IV] that the universal functions for observables are insensitive to the behavior 
of IT at k, oo ~ A. In view of this, we use in Sec |V| and Sec the result for Nq/ p s obtained 
from ( |D19| ) in the case when the polarization operator is computed without a cutoff in the 
momentum and frequency integration. Then U(k,uj) is given by (pil|) and performing the 
integration in (pi8|), we find 
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Ps N 



3tt 2 N 



log 



' NA 
16 p s 



(D20) 



This equation we use in (|4.2|) . 

We will also need the result for p s expressed in terms of the spin-stiffness at iV = oo. 
From flDTJ) ([DTOD and flDI5| ), we obtain 



= P. 



N=oo 



1 + 



4w 2 



N J (2tt)3 n*(Mw) + ^ 2 ) 



;D2i) 



Finally, we deduce from ( p21 ) the critical exponent for p s . For this, we perform the 
integration in (|D21| ) with the logarithmic accuracy using ( pil| ) and (P5|) , and exponentiate 
the result. We then obtain 



Ps 



, 9c ~ 9 y 



(D22) 



where 



2(3 



32 



1+T] 3tt 2 N 
and rj = 8/(37r 2 iV) is the critical exponent for spin correlations at g c . 



(D23) 



74 



REFERENCES 



1 S. Chakravarty in High- Temperature Superconductivity, eds. K. Bedell, D. Coffey, D.E. 
Meltzer, D. Pines and J.R. Schrieffer, Addison- Wesley, p. 136 (1990). 

2 E.Manousakis, Rev. Mod. Phys, 63, 1 (1991). 

3 S. Sachdev, in Low Dimensional Quantum Field Theories for Condensed Matter Physi- 
cists, Proceedings of the Trieste Summer School 1992, World Scientific, to be published, 
and references therein. Available as paper 9303014 on cond-mat@babbage.sissa.it. 

4 A.P. Ramirez, G.P Espinosa, and A.S. Cooper, Phys. Rev. Lett. 64, 2070 (1990). 

5 C. Broholm, G. Aeppli, G.P. Espinosa, and A.S. Cooper, Phys. Rev. Lett. 65, 3173 

(1991) . 

6 G. Aeppli, C. Broholm, and A. Ramirez in Proceedings of the Kagome Workshop, NEC 

Research Institute, Princeton, NJ (unpublished). 
7 V. Elser, Phys. Rev. Lett., 62, 2405 (1990). 

8 S.J. Clarke, A. Harrison, T.E. Mason, G.J. Mclntyre and D. Visser, J. Phys. Cond. 
Matter, 4, L71 (1992). 

9 For a review on recent experiments see e.g. A. Harrison, in "Annual Reports on the 
Progress of Chemistry (Royal Soc. of Chemistry, UK), 87A, 211, (1992); 88 A, 447, 

(1992) . 

10 V.L. Pokrovskii, Adv. Phys. 28, 595 (1979). 

11 A.F. Andreev and V. I. Marchenko, Sov. Phys. Usp. 23, 21, (1980). 

12 S. Chakravarty, B.I. Halperin and D.R. Nelson, Phys. Rev. Lett. 60, 1057 (1988); Phys. 
Rev. B 39, 2344 (1989). 

13 S. Tyc, B.I. Halperin, S. Chakravarty, Phys. Rev. Lett. 62, 835 (1989). 

14 V. Kalmeyer and R.B. Laughlin, Phys. Rev. Lett., 59, 2095 (1987); X.G. Wen, F. Wilczek 
and A. Zee, Phys. REv. B 39, 11413 (1989). 

15 N. Read and S. Sachdev, Phys. Rev. Lett. 62, 1694 (1989); Phys. Rev. B 42, 4568 (1990). 

16 N. Read and S. Sachdev Phys. Rev. Lett. 66, 1773 (1991); S. Sachdev and N. Read, Int. 
J. Mod. Phys. B5, 219 (1991). 

17 P.Chandra, P. Coleman and A.I. Larkin, J. Phys. Cond. Matter 2, 7933 (1990). 

18 I. Affleck and J. Brad Marston, Phys. Rev. B 37, 3774, (1988). 

19 PB. Wiegmann, Phys. Rev. Lett. 60, 821, (1988); D.V. Khveshchenko and PB. Wieg- 
mann, Mod. Phys. Lett. B 4, 17 (1990). 

20 PW. Anderson, Science 235, 1196 (1987). 

21 S. Sachdev and Jinwu Ye, Phys. Rev. Lett. 69, 2411 (1992). 

22 A. V. Chubukov and S. Sachdev, Phys. Rev. Lett. 71, 169 (1993); erratum, to be pub- 
lished. 

23 B. Andraka and A.M. Tsvelik, Phys. Rev. Lett. 67, 2886 (1991). 

24 A.J. Millis, preprint; J.A. Hertz, Phys. Rev. B 14, 525 (1976). 

25 P. Hasenfratz and F. Niedermayer, Phys. Lett. B 268, 231 (1991); University of Bern 
preprint BUTP-92/46. 

26 P. Hasenfratz, M. Maggiore and F. Niedermayer, Phys. Lett. B 245, 522 (1990); P. 
Hasenfratz and F. Niedermayer, Phys. Lett. B 245, 529 (1990). 

27 S. Chakravarty and R. Orbach, Phys. Rev. Lett 64, 224 (1990). 



75 



28 Y. Endoh et. al, Phys. Rev. B 37, 7443 (1988); K. Yamada et. al, Phys. Rev. B 40, 
4557 (1989). 

29 S.M. Hayden, G. Aeppli, H. Mook, D. Rytz, M.F. Hundley, and Z. Fisk, Phys. Rev. Lett. 
66, 821 (1991); S.M. Hayden, G. Aeppli, R. Osborn, A.D. Taylor, T.G. Perring, S.-W. 
Cheong, and Z. Fisk, Phys. Rev. Lett. , 67, 3622 (1991). 

30 T. Imai, CP Slichter, K. Yoshimura and K. Kosuge, Phys. Rev. Lett 70, 10002, (1993). 
31 D.C. Johnson, Phys. Rev. Lett. 62, 957 (1989). 

32 R.R.P. Singh and M. Gelfand, Phys. Rev. B 42, 966 (1990). 

33 H.Q. Ding and M. Makivic, Phys. Rev. Lett, 64, 1449 (1990); Phys. Rev. B 43, 3662 
(1990). 

34 A. Sokol, S. Bacci and E. Gagliano, Phys. Rev. B 47, 14646 (1993). 

35 B. Keimer, R.J. Birgeneau, A. Cassanho, Y. Endoh, R.W. Erwin, M.A. Kastner, and G. 
Shirane, Phys. Rev. Lett. 67, 1930 (1991). 

36 B. Keimer, N. Belk, R.J. Birgeneau, A. Cassanho, C.Y. Chen, M. Greven, M.A. Kastner, 
A. Aharony, Y. Endoh, R.W. Erwin, and G. Shirane, Phys. Rev. B 46, 14034 (1992). 

37 CM. Varma P.B. Littlewood, S. Schmitt-Rink, E. Abrahams, and A.E. Ruckenstein, 
Phys. Rev. Lett. 63, 1996 (1989). 

38 L.B. Ioffe and A.I. Larkin, Int. J. Mod. Phys. B 2, 203 (1988). 
39 A.V. Chubukov, Phys. Rev. B 44, 392 (1991). 

40 B.I. Halperin and W.M. Saslow, Phys. Rev. B 16, 2154 (1977). 

41 T. Dombre and N. Read, Phys. Rev. B 39, 6797 (1989). 

42 P. Azaria, B. Delamotte, and T. Jolicoeur, Phys. Rev. Lett. 64, 3175 (1990). 

43 A.B. Harris, C. Kallin and A.J. Berlinsky, Phys. Rev B 45, 2889 (1992); J.T. Chalker, 
PS. Holdsworth and E.F. Shender, Phys. Rev. Lett. 68, 855 (1992); P. Chandra, P. 
Coleman and I. Ritchey, Phys. Rev. B,to appear; A.V. Chubukov, Phys. Rev. Lett. 69, 
832 (1992); J. von Delft and CL. Henley, Phys. Rev. Lett. 69, 3236 (1992); J.N. Reimers, 
A.J. Berlinsky and A.-C Shi, Phys. Rev. B 43, 865 (1991); R.R.P. Singh and D. Huse, 
Phys. Rev. Lett. 68, 1706 (1992); D.L. Huber and WY. Cheong, Phys. Rev. B 47, 3220 
(1993). 

44 S. Sachdev, Phys. Rev. B 45, 12377 (1992). 

45 S-k. Ma, Modern Theory of Critical Phenomena, Benjamin/Cummings, Reading (1976). 

46 J.B. Parkinson, J. Phys. C 2, 2012 (1969); B.S. Shastry and B. Shraiman, Phys. Rev. 
Lett. 65, 1068 (1990); R.R.P. Singh, Comments Cond. Matt. Phys. 15, 241 (1991). 

47 M.E. Fisher, M.N. Barber, and D. Jasnow, Phys. Rev. A 8, 1111 (1973). 
48 B.D. Josephson, Phys. Lett. 21, 608 (1966). 

49 A.H. Castro Neto and E. Fradkin, paper 9301009 on cond-mat@babbage.sissa.it 

50 A.N. Vasil'ev, Yu.M. Pis'mak and Yu.R. Honkonen, Teor. Mat. Fiz. 46, 157 (1981). 

51 C. Holm and W. Janke, preprint; P. Peczak, A.M. Ferrenberg, and D.P. Landau, Phys. 
Rev. B 43, 6087 (1991). 

52 In an earlier version of the manuscript we had an incorrect conjecture of the form of the 
subleading term in the large x± t 2 expansion. S. Chakravarty pointed out to us, and we 
determined by explicit computation the present correct form. 

53 D. Stauffer, M. Ferer, and M. Wortis, Phys. Rev. Lett. 29, 345 (1972). 

54 P.C Hohenberg, A. Aharony, B.I. Halperin, and E.D. Siggia, Phys. Rev. B 13, 2986 
(1976); C. Bervillier, ibid. 14, 4964 (1976). 



76 



55 M.P.A. Fisher, G. Grinstein, and S.M. Girvin, Phys. Rev. Lett. 64, 587 (1990); K. Kim 
and P.B. Weichmann, Phys. Rev. B 43, 13583 (1991); M.-C. Cha etal, Phys. Rev. B 44, 
6883 (1991). 

56 M.N. Barber in Phase Transitions and Critical Phenomena, ed. C. Domb and J. Lebowitz, 
v. 8, pg. 145, Acad. Press (New York). 

57 V. Privman and M.E. Fisher, Phys. Rev. B 30, 322 (1984). 
58 D.S. Fisher, Phys. Rev. B 39, 11783 (1989). 

59 M.E. Fisher and P.-G. de Gennes, C.R. Acad. Sci. Ser. B 287, 207 (1978). 
60 A.B. Zamalodchikov, Pis'ma Zh. Eksp. Teor. Fiz 43, 565 (1986); [JETP Lett. 43, 730 
(1986). 

61 H.W.J. Blote, J.L. Cardy, and M.P Nightingale, Phys. Rev. Lett. 56, 742 (1986); I. 

Affleck Phys. Rev. Lett. 56, 746 (1986). 
62 D.P. Arovas and D. Auerbach, Phys. Rev. B 38, 316 (1988); Phys. Rev. Lett. 61, 617 

(1988). 

63 B. Rosenstein, B.J. Warr and S.H. Park, Nucl. Phys. B 336, 435 (1990). 
64 E. Brezin, J. de Physique 43, 15 (1982); M. Henkel and C. Hoeger, Z. Phys. B 55, 67 
(1984); S. Singh and R.K. Pathria, Phys. Rev. B 31, 4483 (1985) and references therein. 
65 R.F. Dashen, S-K. Ma and R. Rajaraman, Phys. Rev. D 11, 1499 (1975). 

66 S. Sachdev, Phys. Lett. B 309, 285 (1993). 

67 A.M. Polyakov, Gauge Fields and Strings, Harwood, New York (1987). 

68 S. Coleman, R. Jackiw and D. Politzer, Phys. Rev. D 10, 2491 (1974). 

69 LA. Aref'eva, Ann. Phys. (N.Y.) 117, 393 (1979). 

70 J.L. Cardy, J. Phys. A. 17, L385 (1984); R. Shankar and S. Sachdev, unpublished. 

71 A.M. Polyakov, Phys. Lett. B 59, 79 (1975). 

72 P.B. Wiegmann, Pis'ma Zh. Eksp. Teor. Fiz. 41, 79 (1985) [JETP Lett. 41, 95 (1985). 

73 E. Brezin and J. Zinn- Justin, Phys. Rev. B 14, 3110 (1976). 

74 A.M. Polyakov and P.B. Wiegmann, Phys. Lett. B 131, 121 (1983). 

75 A. V. Chubukov, Phys. Rev. B 44, 12318 (1992). 

76 M.I. Kaganov and A.V. Chubukov, Usp. Fiz. Nauk 153, 537 (1987) [Sov. Phys. Usp. 30, 
1015 (1987)]; in "Spin Waves and Magnetic Dielectrics" eds. A.S. Borovik-Romanov and 
S.K. Sinha, Elsevier Science Publ. (1988). 

77 S. Tyc and B.I. Halperin, Phys. Rev. B 42, 2096 (1990). 

78 B.I. Halperin and PC. Hohenberg, Phys. Rev. 177, 952 (1969). 

79 In the recent paper by one of us^l, it was suggested that the damping may become compa- 
rable to the real part of the quasiparticle energy at the spatial scales which parametrically 
exceed £. The more sophisticated analysis presented here shows that it is more likely that 
2D antiferromagnet has the same typical spatial scale (correlation length) for both static 
and dynamic phenomenon, as is predicted by the dynamical scaling hypothesis. 

80 J. Igarashi, Phys. Rev. B 46, 10763 (1992). 

81 R.R.P. Singh, Phys. Rev. B 39, 9760 (1989); see also R.R.P. Singh and D. Huse, Phys. 
Rev. B 40, 7247 (1989). 

82 A. Singh and Z. Tesanovic, Phys. Rev. B 41, 614 (1990); A. Singh, Z. Tesanovic and J.H. 
Kim, Phys. Rev. B 44, 7757 (1991); A. Auerbach and B.E. Larson, Phys. Rev. B 43, 
7800 (1991). 

83 B.I. Shraiman and E.D. Siggia, Phys. Rev. Lett. 61, 467 (1988); Phys. Rev. B 42, 2485 



77 



(1990). 

84 B.I. Shraiman and E.D. Siggia, Phys. Rev. B 46, 8305 (1992). 

85 A. V. Chubukov and D. Frenkel, Phys. Rev. B 46, 11884 (1992). 

86 J. Gan, N. Andrey and P. Coleman, J. Phys. Condens. Matter 3, 3537 (1991). 

87 S. Sachdev, unpublished. 

88 R.R.P. Singh, Phys. Rev. B 39, 9760 (1989). 

89 A. Millis and H. Monien, Phys. Rev. Lett. 70, 2810 (1993). 

90 D.C. Johnson, S.K. Sinha, A.J. Jacobson and J.M. Newsam, Physica (Amsterdam), 153- 
155 C, 572 (1988). 

91 F. Mila and T.M. Rice, Physica C 157, 561 (1989); S. Shastry Phys. Rev. Lett., 63, 1288 
(1989). 

92 A.J. Millis, H. Monien and D. Pines, Phys. Rev. B 42, 167 (1990). 

93 E. Manousakis and R. Salvador, Phys. Rev. B 39, 575 (1989); E. Manousakis, Phys. Rev. 
B 45, 7570, (1992). 

94 E. Manousakis and R. Salvador, Phys. Rev B 40, 2205 (1989). 

95 A. Sokol, private communication. 

96 A. Chubukov, S. Sachdev and A. Sokol, Phys. Rev. B, submitted. 
97 M. Makivic and M. Jarrell, Phys. Rev. Lett., 68, 1770 (1992). 

98 M.P.A. Fisher, P.B. Weichmann, G. Grinstein, and D.S. Fisher, Phys. Rev. B 40, 546 
(1989). 

"R.N. Bhatt and PA. Lee, Phys. Rev. Lett. 48, 344 (1982). 

100 I. Affleck and A.W.W. Ludwig, Nucl. Phys. B360, 641 (1991). 

101 K. Binder and A.P Young, Rev. Mod. Phys. 58, 801 (1986). 

102 G. Murthy and S. Sachdev, Nucl. Phys. B344, 557 (1990). 

103 D.S. Fisher, private communication. 

104 U. Wolff, Phys. Rev. Lett. 62, 361 (1989). 



78 



FIGURES 

FIG. 1. Phase diagram of TL (Eqn. fll-lQ ) as a function of g and temperature T (after 
Refit). The coupling g measures the strength of the quantum fluctuations. It is inversely pro- 
portional to S for large spin and its value also depends on the ratios of the Jij. The parameters 
xi = NkBT/(2irp s ), and x<i = ksT/A control the scaling properties of the antiferromagnet (here 
p s is the spin-stiffness of the Neel-ordered ground state, and A is the spin 1 gap in the quan- 
tum-disordered ground state). 

FIG. 2. Properties of the nearly-critical antiferromagnet as a function of the observation 
wavevector k, or frequency u> in the three different regions of Fig. |l[ The appropriate regime 
is determined by the larger of tick/ {ksT) or huj/ksT. In the renormalized-classical regime, £ is 
the actual correlation length, while £j is a Josephson correlation length related to the spin-stiffness 
by hc/^j = p s /T, with T a universal number. In the quantum disordered region, A is the gap for 
S = 1 excitations at T = 0. The thermodynamic behavior in the various regions is discussed in 
the text. 

FIG. 3. The scaling function Im$i s (A;,a;, oo) for the staggered susceptibility in the quan- 
tum-critical region. The results have been computed in a 1/N expansion to order 1/N and evaluated 
for N = 3. The shoulder on the peaks is due to a threshold to three spin-wave decay. 



FIG. 4. Scaling function Hi(fc,oo) for the structure factor in the quantum-critical region. 

FIG. 5. Scaling function Fi(p, oo) for the imaginary part of the local susceptibility in the 
quantum-critical region. The oscillations at large ZJ are due to a finite step-size in the momentum 
integration. 

FIG. 6. Feynman graph for the staggered and uniform susceptibilities to order 1/N. The solid 
line represents the ri£ field and the dashed line represents the propagator of the A field (polarization 
operator), where A is the Lagrange multiplier imposing the constraint. 

FIG. 7. Quantum Monte-CarkU (squares) and our theoretical (line) results for the uniform 
susceptibility x u = (3J(ah/gp,B) 2 )Xu °f a square lattice spin-1/2 Heisenberg antiferromagnet (a 
is the lattice spacing). The experimental results for weakly doped La2CuO^ are very close to the 
Monte-Carlo data£x There are no adjustable parameters in the theoretical result ( |1.24| ). Over 
the range of T plotted, the function Qi(xi) (recall x\ = NkBT/(2ivp s ) ) is very close to its large 
x\ behavior given in Eqn. ( p,35 , 4~T3| ). We used the theoretical result at N = 3. The theoretical 



and experimental slopes agree remarkably well. The good agreement in the intercept is somewhat 
surprising as its theoretical value is known only at N = oo. 
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FIG. 8. Theoretical scaling function at N = oo, Eqn. (|7.6|) (line) and Quantum Monte-Carlo 
datall (squares) for the equal-time structure factor in the 5 = 1/2 Heisenberg antiferromagnet at 
ksT = 0.42J. The correlation length is taken from the Monte-Carlo data (£ ~ 5.7a). We expect 
that for most of the values of k£ plotted, the antiferromagnet is in the quantum-critical regime 
where 1/N corrections to (|7.6| ) are small. The only adjustable parameter in the theoretical curve is 
the A:— independent overall factor S in (7.6). The best fit value of S was found to be ~ 1.45 times 
larger than our N = oo result. 



FIG. 9. Schematic of the renormalization group flows of the action S (Eqn. (|B2|)) of a classical 
XY model on a cubic lattice. The coupling h& is a cubic anisotropy perturbation which is irrelevant 
at the critical fixed point T = T c , /14 = 0. Its neglect is however dangerous in the low temperature 
phase because is relevant at the T = 0, /14 = fixed point. 
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TABLES 



TABLE I. Results of the Monte Carlo simulation of the classical statistical-mechanics model 
(C2) at K = K c = 0.6930. We used a box of L x L x L T sites with periodic boundary conditions. 
The stiffness p T is defined in (|C4|) . The last column is obtained by a polynomial extrapolation 
in inverse powers of 1/L. The three runs had 70000, 70000, 210000 flips per spin respectively. A 
weighted average of the three runs was used in the extrapolation. 





L 




L T Pr 




liniL^oo L T p T 






Run 1 


Run 2 


Run 3 






10 


0.3983 


0.3866 


0.3898 






15 


0.3693 


0.3701 


0.3663 




5 


20 


0.3596 


0.3515 


0.3527 


0.3257 




25 


0.3511 


0.3574 


0.3478 






30 


0.3483 


0.3529 


0.3511 






10 


0.4138 


0.4087 


0.4099 






15 


0.3718 


0.3766 


0.3747 




7 


20 


0.3650 


0.3625 


0.3506 


0.3037 




25 


0.3529 


0.3419 


0.3442 






30 


0.3424 


0.3405 


0.3402 






15 


0.4079 


0.4049 


0.4020 




10 


20 


0.3861 


0.3743 


0.3713 


0.2890 




25 


0.3650 


0.3671 


0.3511 






30 


0.3548 


0.3521 


0.3425 
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